Singlet  Delta  Oxygen:  A  Quantitative  Analysis  Using 
Off-Axis  Integrated-Cavity-Output-Spectroscopy  (ICOS) 


THESIS 

Jeffrey  E  Gallagher,  Captain,  USAF 
AFIT  /  GAP /ENP /06-05 


DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 

Wright-Patterson  Air  Force  Base,  Ohio 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
policy  or  position  of  the  United  States  Air  Force,  the  Department  of  Defense  or  the 
United  States  Government. 


AFIT/GAP /ENP/06-05 


Singlet  Delta  Oxygen:  A  Quantitative  Analysis  Using 
Off-Axis  Integrated-Cavity-Output-Spectroscopy  (ICOS) 


thesis 


Presented  to  the  Faculty 
Department  of  Engineering  Physics 
Graduate  School  of  Engineering  and  Management 
Air  Force  Institute  of  Technology 
Air  University 

Air  Education  and  Training  Command 
In  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Master  of  Science  (Applied  Physics) 


Jeffrey  E  Gallagher 
Captain,  USAF 


March  2006 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


AFIT/GAP /ENP/06-05 


Singlet  Delta  Oxygen:  A  Quantitative  Analysis  Using 
Off-Axis  Integrated-Cavity-Output-Spectroscopy  (ICOS) 


Jeffrey  E  Gallagher 
Captain,  USAF 


Approved: 


/signed/  03  Mar  2006 

Glen  P.  Perrarn  (Chairman)  date 

/signed/  03  Mar  2006 

Paul  J.  Wolf  (Member)  date 

/signed/  03  Mar  2006 

David  E.  Weeks  (Member)  date 

/signed/  03  Mar  2006 


Skip  Williams  (Member) 


date 


AFIT /GAP /ENP/06-05 


Abstract 


A  new  spectroscopic  technique  applicable  to  the  detection  of  ultra-weak  and  for¬ 
bidden  molecular  transitions  is  presented.  The  method  is  based  on  off-axis  integrated- 
cavity-output  spectroscopy  (ICOS)  and  has  been  applied  to  the  detection  of  singlet 
(a1A9)  oxygen  in  the  (1,0)  band  of  the  clcctric-quadrupolc  allowed  Noxon  system 
( b 1E+  <—  a1As)  of  oxygen.  The  details  of  the  method  as  well  as  spectroscopic  data 
confirming  the  absolute  line  positions  of  twenty-one  lines  in  the  (1,0)  band  of  the 
a1  A „  have  been  completed  to  within  0.004  cm"1.  A  Boltzmann  analysis 


6'y+ 


verified  previously  calculated  and  unproven  integrated  cross  sections  for  each  of  the 
lines.  Ten  lines  were  also  pressure  broadened  from  25  Torr  to  100  Torr,  and  pres¬ 
sure  broadening  coefficients  were  determined  from  =  3.61  x  10~6  cm"1  /torr  to 
5.81  x  10~6  cm"1  /torr.  These  coefficients  illustrate  that  rotationally-inelastic  colli¬ 
sions  are  the  dominant  mechanism  in  molecular  pressure  broadening.  The  singlet  oxy¬ 
gen  was  generated  in  a  microwave  plasma  and  the  afterglow  passed  through  an  off-axis 
ICOS  measurement  system  consisting  of  an  82  cm  long,  high-finesse  ( F  ss  248000) 
optical  cavity.  The  mirror  reflectivity,  R,  was  determined  by  performing  a  cavity- 
ringdown  measurement  and  observing  ringdown  times  of  220  /zs  in  a  range  from  1494 
nm  to  1512  nm.  The  diode  laser  was  current  tuned  over  ca.  0.33  cm-1  at  5-25 
Hz.  Light  exiting  the  cavity  was  focused  onto  an  InGaAs  detector,  and  the  cavity 
transmission  was  recorded  as  a  function  of  laser  frequency.  The  absorption  within 
the  cavity  leads  to  a  decrease  in  transmitted  intensity,  A/,  given  by  =  GA 


Io  l+GA  ’ 

where  ID  is  incident  laser  intensity,  A  is  the  single-pass  absorption  (f«  5  x  10  8),  and 


G  = 


2  R2 

1  -R2 


80,  000  is  the  cavity  enhancement  factor. 
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Singlet  Delta  Oxygen:  A  Quantitative  Analysis  Using 
Off-Axis  Integrated-Cavity-Output-Spectroscopy  (ICOS) 


I.  Introduction 

The  61Es+ "  “'As  system  of  oxyEen  is  an  electric_quadrupole  transition  md  is 

JL  therefore  very  weak  and  extremely  difficult  to  detect.  In  addition,  this  partic¬ 
ular  absorption  band  (referred  to  as  the  Noxon  Band  after  John  Noxon  [15])  has  a 
radiative  lifetime  of  about  72  minutes  [7] ,  which  classifies  this  as  a  metastable  species. 
Measuring  the  absorption  of  a  metastable  species  is  in  itself  very  difficult.  This  weak 
signal  combined  with  the  long  lifetime  is  why  these  transitions  have  not  been  directly 
observed  until  quite  recently  [24],  The  research  goals  for  this  experiment  and  the 
benefits  to  the  Department  of  Defense  (DoD)  are  explained  in  the  following  sections. 

1.1  Research  Goals 

The  primary  goal  for  this  research  effort  is  to  utilize  a  relatively  new  spectro¬ 
scopic  technique  and  demonstrate  it’s  ability  to  provide  quantitative  data  of  singlet 
delta  oxygen.  This  thesis  will  focus  on  three  areas  of  characterization  to  achieve  this 
goal.  First,  the  absolute  line  positions  will  be  determined  and  compared  to  values 
derived  from  the  most  recent  theory.  Second,  the  integrated  absorption  cross-sections 
will  be  verified  using  Boltzmann  analysis.  Finally,  pressure  broadening  coefficients 
will  be  determined  for  select  lines  using  a  range  of  pressures  from  25  to  100  Torr. 
After  the  coefficients  are  determined,  this  experiment  will  hope  to  be  the  first  in  dis¬ 
cerning  a  rotational  dependence  on  pressure  broadening  for  an  electric-quadrupole 
allowed  transition. 

1.2  Benefits  to  the  DoD 

The  DoD  is  concerned  about  characterizing  singlet  delta  oxygen  because  this 
knowledge  will  help  diagnose  the  collisional  environment  of  the  Chemical  Oxygen- 
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Iodine  Laser  (COIL).  Additional  areas  that  singlet  delta  oxygen  (02(a))  plays  a  role 
in  include  biological  weapons  detection,  satellite  communications,  semiconductor  tech¬ 
nology,  and  propulsion  systems. 

1.2.1  The  Chemical  Oxy gen- Iodine  Laser.  The  COIL  as  a  mega- watt  class 
chemical  laser  that  is  currently  being  integrated  onto  an  aircraft  for  the  purpose  of 
missile  defense  called  the  Airborne  Laser  program  or  ABL.  The  production  of  singlet 
delta  oxygen  [02(a1A3)],  or  simply  02(a),  is  a  very  important  part  of  the  lasing 
process  in  the  COIL  since  this  is  the  molecule  used  as  the  energy  reservoir  due  to 
it’s  long  lifetime.  As  the  02(a)  goes  through  the  supersonic  nozzle,  gaseous  Iodine 
is  introduced  into  the  system  where  the  02(a)  eventually  transfers  it’s  energy  in  a 
near-resonant  collision  to  a  dissociated  Iodine  atom  (J(2P3/2)).  The  actual  lasing 
transition  occurs  between  /(2Pi/2)  and  /(2P3/2),  in  the  following  process, 

J(2P1/2)  +  h  u  — *  J(2P3/2)  +  2  hv  (1.1) 

02(a),  however,  plays  a  crucial  role  in  a  series  of  chemical  reactions  which  lead  to  the 
production  of  the  needed  excited  Iodine  atom,  /(2Pi/2)  [13]. 

n02(a)  +  /2  — »  /(2P3/2)  +  J(2P3/2)  +  n02(X3T,g)  (1-2) 

02(a)  +  J(2P3/2)  J(2P1/2)  +  02(X%)  (1.3) 

The  n  is  used  in  equation  1.2  because  more  than  one  02(a)  molecules  are  required 
to  disassociate  the  Iodine  molecule.  In  equation  1.3  the  02(a)  molecule  reacts  with 
I(2P 3/2)  in  a  near  resonant  transfer  to  produce  ground  state  oxygen  02(A")  and  the 
desired  /(2Pi/2).  It  is  therefore  very  important  to  understand  the  properties  of  02(a) 
in  order  to  better  control  its  production.  Scientists  and  engineers  in  the  COIL  com¬ 
munity  are  striving  for  a  method  to  directly  measure  the  concentration  of  02(a)  prior 
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to  the  lasing  process.  The  singlet  oxygen  generator  (SOG)  used  in  the  current  COIL 
is  believed  to  produce  high  yields  of  02(a )  molecules  because  of  the  following  reaction. 
It  produces  only  02(a).  The  off-axis  ICOS  apparatus  used  in  this  experiment  will  be 
able  to  quantify  the  number  density  of  this  02(a)  production  prior  to  entering  the 
nozzle  assembly.  The  research  and  technique  used  in  this  thesis  will  help  develop  that 
method. 

1.2.2  Atmospheric  Modeling.  In  the  upper  atmosphere,  namely  the  ther¬ 
mosphere,  singlet  delta  oxygen  is  fairly  abundant  and  contributes  strongly  to  the 
radiation  transfer  which  takes  place  there.  Much  like  the  COIL,  the  kinetics  involved 
in  the  thermosphere  are  not  completely  understood.  Consequently,  gaining  a  better 
understanding  of  the  earth’s  atmosphere  will  be  very  important  and  would  improve 
our  satellite  communications  network  and  other  forms  of  energy  propagating  through 
the  atmosphere. 

1.2.3  Other  Benefits.  O11  the  topic  of  biological  weapons  detection,  it  has 
been  shown  that  air  plasmas  produce  large  quantities  of  02(a)  with  demonstrated 
reactivity  with  hydrocarbons  and  biological  agents.  In  the  semiconductor  industry, 
02(a)  and  atomic  oxygen  are  used  for  ashing  photoresist  films  onto  semiconductor 
wafers.  Finally,  the  study  of  02(a)  may  also  lead  to  more  efficient  propulsion  systems. 
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II.  Background 


Absorption  from  02(a)  to  singlet  sigma  oxygen  (02(b))  is  very  weak  and  the  op¬ 
tical  transitions  are  highly  forbidden.  This  chapter  will  explore  the  underlying 
reasons  why  this  molecule  is  difficult  to  detect  and  quantify.  First,  a  short  history 
will  be  presented  of  those  who  tried  to  characterize  02(a1Ag )  both  experimentally  and 
theoretically.  Then,  various  topics  on  the  theory  will  be  discussed.  These  topics  are 
intended  to  aid  the  reader  in  comprehending  the  data  contained  in  this  thesis.  Finally 
there  is  a  short  discussion  on  the  use  of  the  Hamiltonian  to  derive  the  theoretical  line 
positions. 

2.1  History  of  Detecting  the  <—  a1Ag  Transition  of  Oxygen 

The  history  of  detecting  the  <—  alAg  system  of  molecular  oxygen  is  very 

brief.  It  began  when  Noxon  [15]  first  observed  the  b — ►  a1Ag  emission  spectra 
in  a  laboratory  in  1961.  Constrained  with  1950s  technology,  Noxon  merely  saw  a 


Figure  2.1:  John  Noxon’s  emission  spectra  of  the  transition  from  02{b)  to  02(a).  [15] 

convolution  of  the  Q  branch  viewed  as  a  single  peak.  None  of  the  rotational  levels 
were  discernable,  although  he  noted  their  theoretical  positions  on  his  spectra.  With 
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the  inability  to  achieve  better  resolution  of  these  lines,  the  theoreticians  then  played  a 
part  to  better  understand  these  highly  forbidden  systems.  In  1964,  Ying-Nan  Chiu  [3] 
derived  the  theoretical  line-strength  formulas  for  both  magnetic-dipole  and  electric- 
quadrupole  radiation.  Then,  in  1986,  25  years  after  John  Noxon’s  work,  German 
and  Canadian  scientists  [6]  unearthed  Noxon’s  and  Chiu’s  work  and  calculated  and 
detected  twenty  rotational  levels  in  each  of  the  five  allowed  branches  of  the  emission 
of  the  b 1S+  <—  a1  A g  transition.  Fink  et  al  was  the  first  to  validate  the  relative  signal 


Figure  2.2:  Group  from  Germany  and  Canada  detected  all  five  branches  using  a 

fourier  transform  spectrometer.  They  achieved  a  resolution  of  0.2cm”1  and  were  the 
first  to  validate  the  relative  signal  strength  intensities  among  the  rotational  lines.  [6] 


strength  intensities  among  the  rotational  levels  in  each  of  the  five  branches. 

Another  15  years  go  by  before  anyone  else  detects  this  transition.  A  group 
from  Russia  [18]  in  2001  detected  the  absorption  spectra  of  the  Noxon  system  for 
the  first  time  with  a  technique  called  intra-cavity  laser  spectroscopy.  Finally  in  2003, 
Williams  et  al  further  enhanced  the  off-axis  ICOS  method  and  quantitatively  detected 
rotational  levels  in  the  Q  branch  of  the  b1^  <—  a1Ag  system  for  the  (0,0)  band  [13] 
and  the  (1,0)  band  [24], 
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Figure  2.3:  The  Russians  were  the  first  to  detect  the  absorption  lines  of  the  Noxon 
system.  They  used  a  technique  called  intra-cavity  absorption  spectroscopy  and  they 
achieved  a  resolution  of  0.018cm_1.  [18] 
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Figure  2.4:  Williams  et  al  was  the  first  to  detect  the  (1,0)  band  of  the  Q  branch. 
Williams  used  O’Keefe’s  off-axis  1COS  technique  to  finely  resolve  the  rotational 
lines.  [24] 
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Figure  2.5:  Williams  et  al  used  the  off-axis  ICOS  technique  to  finely  resolve  the 
(0,0)  band  of  the  Q  branch.  [13] 
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2.2  History  of  Cavity  Ringdown 

In  the  early  90’s  the  cavity  ringdown  technique  emerged  from  the  spectroscopic 
community  [11].  This  technique  allowed  the  experimentalists  to  observe  weaker  tran¬ 
sitions  in  the  spectra  and  set  the  stage  for  the  eventual  detection  and  quantification 
of  singlet  oxygen  concentrations.  In  1992  Romanini  and  Lehmann  [23]  describe  their 
successes  in  detecting  weak  overtones  of  HCN  using  a  ringdown  cavity  (RDC)  absorp¬ 
tion  scheme.  Here  they  achieved  ringdown  times  of  40  ys  and  effective  path  lengths  of 
24  km.  This  technique  became  the  spring  board  for  detecting  very  weak  transitions. 

By  1998,  the  RDC  absorption  scheme,  now  called  cavity  ringdown  spectroscopy 
(CRDS)  had  become  widely  accepted  by  spectroscopists.  In  this  year,  an  extensive 
survey  on  the  theory  and  applications  of  CRDS  was  written  by  scientists  from  the 
UK  [14].  Three  years  later,  radiative  lifetimes  of  02(a1A3)  were  measured  using 
CRDS  [7], 

Although  CRDS  is  not  used  in  this  experiment,  the  technique  in  this  effort 
will  use  the  method  of  obtaining  the  ringdown  time  from  the  optical  cavity.  This 
ringdown  time  will  be  a  key  factor  in  amplifying  the  weak  signals  associated  with  the 
b1^  <—  a1Ag  transition  of  oxygen. 

2.3  History  of  Off-axis  ICOS 

In  1998,  O’Keefe  [16]  explored  a  new  variation  of  CRDS,  called  integrated- 
cavity-output-spectroscopy  (ICOS).  ICOS  employs  the  same  methodology  from  CRDS 
[i.e.  highly  reflective  mirrors,  ringdown  time  measurements  (r)  and  effective  path 
length  measurements,  L],  however,  in  ICOS  instead  of  observing  the  single-pass  out¬ 
put  intensity  (as  is  done  in  CRDS),  all  of  the  output  intensities  are  collected  and 
focused  onto  an  InGaAs  detector  where  the  signals  are  then  integrated  to  produce 
the  absorption  feature  (See  Figure  2.6).  Also,  since  the  effective  path  lengths  are  so 
large  (~  100 km),  Beer’s  Law  can  then  be  represented  as  an  infinite  series.  The  limit 
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at  which  this  series  converges  will  produce  a  factor  that  amplifies  the  weak  signal. 
This  concept  will  be  more  rigorously  explained  in  the  succeeding  chapters. 


Off-axis 

injection 


Figure  2.6:  Here  is  a  basic  schematic  of  the  off-axis  ICOS  technique  .  The  cavity 
outputs  are  collected  and  integrated  to  acquire  the  absorption  feature. 


An  Off-axis  version  of  ICOS  was  then  developed  by  O’Keefe  in  2001.  In  the 
same  year,  Paul  et  al  used  diode  lasers  (instead  of  the  typical  dye  laser)  with  the 
newly  developed  off-axis  ICOS  method  to  get  an  extremely  narrow  bandwidth  [12]. 
A  narrower  bandwidth  means  more  resolving  power,  which  in  turn  allows  for  the 
searching  of  very  weak  absorption  lines.  Paul  et  al  used  the  off-axis  configuration 
because  diode  lasers  are  very  sensitive  to  optical  feedback.  To  rectify  this,  the  diode 
laser  is  simply  injected  off  the  optical  axis.  Not  only  did  the  noise  due  to  the  optical 
feedback  decrease  significantly,  but  the  effective  path  length  increased  by  creating  a 
Lissajous  spot  pattern  on  the  mirrors.  Measurements  in  the  near  infrared  region  using 
off-axis  ICOS  was  then  achievable  [5].  This  research  group  was  also  the  first  to  define 
the  following  relation  (2.1)  which  is  paramount  in  ICOS  theory  today  [12]. 


Al  _  GA 
I(aL  =  0)  “  1  +  GA 


(2.1) 


Where  G  =  ^  12i^2  j  is  the  cavity  enhancement  factor,  A  =  1  —  e~aL  is  the 
single-pass  absorption,  I(aL  =  0)  is  the  laser  intensity  through  the  empty  cavity, 
and  A /  =  I(aL  =  0)  —  I(aL )  is  the  the  change  in  the  steady-state  cavity  output. 
This  equation  is  further  explained  in  Section  3.1.2,  and  Appendix  A  explores  the  full 
mathematical  development  of  this  equation  from  Beers  Law. 


2-4  Physical  Theory 

Singlet  delta  oxygen  is  a  diatomic  homonuclear  molecule  that  couples  its  angular 
momenta  according  to  Hund’s  case  (b).  This  section  will  describe  the  b <—  a1  Ag 
transition,  how  angular  momenta  are  defined  in  molecules,  and  finally  the  symmetry 
properties  and  selection  rules  pertaining  to  the  b 1 E+  a1  Ag  transition. 


2-4-1  The  Transition.  As  mentioned  earlier,  the  Noxon  band  is  the  clectric- 
quadrupole  transition  between  the  first  two  excited  electronic  states  of  molecular 
oxygen  (Figure  2.7).  Specifically,  this  experiment  will  detect  transitions  between  the 
first  vibrational  level  (v  =  0)  of  O 2(a)  and  the  second  vibrational  level  ( v  =  1)  of 
02(b).  This  band  is  written  in  shorthand  as  the  (1,0)  band,  where  the  numbers  refer 
to  the  vibrational  levels  of  the  upper  state  and  lower  state  respectively. 


eV 

8 


6 


4 


2 


0 


Figure  2.7:  The  first  five  electronic  states  of  oxygen.  [2] 
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A  general  statement  to  put  this  transition  into  perspective,  is  the  fact  that  the 
electric-dipole  moment  is  107  times  more  likely  to  occur  than  the  electric-quadrupole 
moment  [8]. 


2.4.2  The  Electric- Quadrupole  Moment.  When  solving  for  the  matrix  ele¬ 
ments  of  the  time-dependent  Schrodinger  equation  for  this  particular  transition,  one 
discovers  that  the  first  non-zero  moment  is  the  electric-quadrupole  moment.  That  is  to 
say,  the  electric-dipole  and  magnetic-dipole  moments  both  vanish  for  this  transition, 
making  the  Noxon  system  a  pure  electric-quadrupole  allowed  transition.  Therefore, 
it  will  be  worthwhile  to  discuss  very  briefly  about  the  electrostatic  field  that  results. 
This  field  is  mathematically  described  by  the  third  term  in  the  multipole  expansion 
of  an  arbitrary  scalar  potential,  0(r),  shown  below  [28], 


N 


N 


<f>(r)  = 


47re0r 


£®+  4^'Eqir‘cose'+ 


2=1 


2=1 


A-Ke^r^ 


N  9 

£2£(3C»>9,-1)  +  ... 


2=1 


(2.2) 


Where  eo  is  the  permittivity  of  free  space,  r  is  the  distance  from  the  origin  to  the 
test  charge,  ry  is  the  distance  from  the  origin  to  some  internal  charge  within  the 
molecule,  qt  is  the  magnitude  of  that  charge  (since  we  are  dealing  with  electrons 
qi  =  6j),  and  6i  is  the  angle  separating  r  and  ry.  It  is  apparent  from  equation  2.2 
that  the  electric-quadrupole  field  weakly  interacts  with  any  test  charge  compared  to 
the  first  two  terms  because  of  it’s  -3  dependence.  After  substituting  the  cos26i  term 
with  direction  cosines,  the  third  summation  in  equation  2.2  becomes  a  tensor.  The 
components  of  the  quadrupole  moment  tensor  can  then  be  written  as  follows  [28], 

N 

Qjk  =  ^2  Qi(3ji  h  -  r2  Sjk)  (2.3) 

2=1 

Where  Qjk  is  the  component  of  the  electric  quadrupole  moment  tensor,  both  jt  and 
hi  spatial  indices,  and  Sjk  is  the  Kronecker  delta  (recall  (j,  k  =  x,y,z )).  It  is  this 
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tensor  which  is  inserted  into  the  three-dimensional  wave  equation  as  an  operator,  Q, 
to  generate  the  scalar  electric-quadrupole  moment  which  is  then  used  to  calculate  the 
integrated  absorption  cross  sections  (oq)  and  the  signal  strength  intensities  (Sj>j"). 
The  quadrupole  term  of  equation  2.2  can  be  written  in  terms  of  the  quadrupole 
moment  tensor  [28], 

<Mr)Q  =  4  r3  2  X]  ^2  IjhQjk  (2.4) 

j=x,y,z  k=x,y,z 

Where  lj  and  4  are  the  direction  cosines.  Since  the  quadrupole  tensor  is  a  symmetric 
tensor1  the  number  of  independent  quantities  reduces  from  nine  to  six.  By  summing 
up  the  diagonal, 


Qxx  +  Qyy  +  Qzz  ~  0  (2-5) 

the  tensor  further  reduces  to  five  terms.  The  tensor  can  be  further  reduced  because 
molecular  oxygen  has  an  axis  of  rotational  symmetry.  This  symmetry  diagonalizes 
the  tensor  because  all  of  the  off  diagonal  terms  cancel  each  other  out,  0  for 

j  ^  k.  Therefore,  only  three  terms  are  left  [28], 

Qxx  +  Qyy  +  Qzz  =  0  (2-6) 

In  fact,  this  symmetry  states  that  for  every  charge  at  a  given  value  x,  there  will  be  an 
equal  charge  at  the  same  numerical  value  for  y  for  the  same  value  of  r.  This  means 
that, 


Ql. 


Q 


yy 


(2.7) 


Qjk  —  Qkj  for  J ^ 
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When  this  relation  is  inserted  into  equation  2.6,  the  result  is  [28] 


Q 


a 

xx 


(2.8) 


Equation  2.8  shows  that  there  is  only  one  independent  component  of  the  quadrupole 
moment  tensor.  This  component  is  typically  written  as  Qa  and  is  defined  as  the  term 
associated  with  the  axis  which  coincides  with  the  axis  of  symmetry.  Therefore,  if  the 
internuclear  axis  of  O2  is  defined  as  the  z-axis,  then  Qa  =  Qazz. 

With  this  definition  of  the  quadrupole  moment,  the  scalar  potential  that  results 
from  an  electric-quadrupole  moment  can  be  written  as  [28], 


0(r)g  = 


Qa  (3 cos29  —  1) 
47reo  4r3 


(2.9) 


Include  some  AP  coupling  references  from  Steinfeld  example  (pg21) 


2-4-3  Angular  Momentum  in  Molecules.  In  defining  the  angular  momentum 
vectors,  diatomic  molecules  must  reference  their  internuclear  axis.  (Figure  2.9  graph¬ 
ically  shows  this  concept.)  The  various  Hund’s  cases  (a-e)  describe  the  ways  in  which 
the  angular  momentum  vectors  are  coupled  to  the  internuclear  axis.  Only  case  (a) 
and  case  (b)  will  be  discussed  here  as  they  relate  to  this  study. 

2-4-3. 1  Hund’s  Case  (a).  The  theoretical  line  positions  in  this  exper¬ 
iment  were  derived  using  a  Hamiltonian  from  Hund’s  case  (a)  for  the  <—  A"3E~ 
and  the  a1  Ag  <—  A3E“  transitions.  The  Hamiltonian  and  spectroscopic  constants 
were  taken  from  a  paper  written  by  Cheah  and  Lee  [27].  They  were  able  to  use  the 
case  (a)  Hamiltonian  because  the  transitions  they  were  looking  at  both  consisted  of 
the  triplet  ground  state  of  molecular  oxygen2. 

Hund’s  case  (a)  states  that  the  electronic  motion  of  the  molecule  (which  includes 
both  the  orbital  and  spin  angular  momentum)  is  strongly  coupled  to  the  internuclear 

2The  term  triplet  refers  to  the  3  degenerate  spin  states 
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axis,  thereby  strongly  defining  the  f2  vector,  which  is  the  total  electronic  orbital  angu¬ 
lar  momentum  about  the  internuclear  axis.  Recall  that  the  total  angular  momentum 
in  atomic  systems  is  defined  as3  J=L+S,  whereas  for  molecules  in  case  (a)  configu¬ 
rations  it  is  written  as  Q  =A+E.  Where  A  is  the  internuclear  axis  component  of  L 
and  E  is  the  internuclear  axis  component  of  S.  Furthermore,  the  Q  and  R  (angular 
momentum  of  nuclear  rotation)  form  the  resultant  J  for  molecules  in  case  (a).  This 
can  all  be  seen  in  Figure  2.8.  Since  this  coupling  configuration  is  very  similar  to  the 


Figure  2.8:  Vector  diagram  for  a  homonuclear  diatomic  [Hund’s  case(a)].  [25] 

symmetric  top,  the  resulting  rotational  energy  levels  for  a  molecule  in  Hund’s  case 
(a)  will  have  the  form  of  the  symmetric  top.  If  we  consider  a  rigid  rotor,  this  energy 
is  [8], 


F(J)  =  BVJ(J  +  1)  +  (C  —  Bv)  A2  (2.10) 

Where  F(J)  is  the  rotational  energy,  Bv  is  the  rotational  constant,  and  C A2  can  be 
considered  as  the  spin-orbit  coupling  term  [25],  but  since  S=0  for  singlet  states,  C=0. 

3Where  L  is  the  electronic  orbital  angular  momentum  and  S  is  the  electronic  spin  angular 
momentum. 
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Therefore  we  have 


F(J)  =  BV[J(J  +  1)  —  A2]  (2.11) 

However,  a  much  better  approximation  is  the  non-rigid  rotor  model,  so  the  equation 
then  becomes, 


F(J)  =  BV[J(J  +  1)  -  A2]  -  DV[J(J  +  1)  -  A2]2  (2.12) 

Where  Dv  is  the  centrifugal  distortion  constant.  Finally,  when  the  vibrational  term 
is  inserted, 


F(J)  =  i/0  +  B0[J( J  +  1)  -  A2]  -  D0[J (J  +  1)  -  A2]2  (2.13) 

This  is  the  final  form  of  the  Hamiltonian  used  by  Cheah  and  Lee  [27] .  They  reasoned 
that  although  the  two  bands  they  looked  at  were  actually  in  transition  from  case  (b) 
to  case  (a),  the  case  (a)  Hamiltonian  gave  the  better  fit  to  their  data.  Therefore, 
in  this  experiment  this  same  Hamiltonian  was  used  to  calculate  the  band  energies 
for  the  <—  a1  Ag  system.  This  was  done  by  simply  subtracting  the  two  term 

energies  found  by  the  case(a)  Hamiltonian.  A  MatLab®  program  was  written  for  this 
experiment  to  generate  a  list  of  band  energies  for  FY,+  <—  a1  A g  (1,0).  This  program 
is  presented  in  Appendix  B. 

2-4-3. 2  Hund’s  case  (b).  The  Hund’s  case  (b)  Hamiltonian  more 
accurately  predicts  the  theoretical  energies  for  the  b 1£+  <—  a1Ag  band,  because  both 
are  singlet  states.  A  singlet  state  has  no  spin  component  and  consequently  f2  is  no 
longer  defined.  Case  (b)  applies  to  all  transitions  where  the  electron  spin  angular 
momentum  S  is  either  weakly  coupled  to  the  internuclear  axis,  or  is  not  present  at 
all  (Refer  to  Figure  2.9).  Therefore,  the  good  quantum  number  is  K,  which  is  the 
total  angular  momentum  without  spin.  The  resulting  equation  for  the  total  angular 
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momentum  is  J  =K+S.  For  the  case  of  02(a)  and  02(b) ,  J  =  K  since  S  =  0.  This 
fact  is  important  when  considering  selection  rules  (See  Figure  2.10). 


Figure  2.9:  Vector  diagram  for  a  homonuclear  diatomic  [Hund’s  case(b)].  [25] 

The  case  (b)  Hamiltonian  was  not  used  in  this  study  because  the  spectroscopic 
constants  (such  as  Bv,  Dv ,  oje\e  and  a jeye)  are  not  yet  available  for  the  b <—  a1  Ag 
transition.  However,  Williams  did  utilize  the  case  (b)  Hamiltonian  in  previous  work 
to  determine  the  line  strength  intensities  (Sj'j")  and  the  integrated  absorption  cross 
sections  (a).  These  calculated  cross  sections  were  used  in  this  study. 

Furthermore,  the  atomic  notation  2S+1L  is  now  written  as  2S+1A  for  molecules. 
Therefore,  singlet  delta  oxygen  (a1  A g)  refers  to  S  =  0  and  A  =  2  and  singlet  sigma 
oxygen  ( b 1E+)  refers  to  S  =  0  and  A  =  0. 

2-4-4  Symmetry  Operations  and  Properties.  The  02(a )  molecule  is  both 
diatomic  and  homonuclear,  therefore  it  is  classified  as  a  Dock  point  group  [2],  The 
Booh  group  is  unique  because  in  addition  to  the  usual  rotation  a  and  reflection  Cn 
operations  for  diatomic  molecules,  the  inversion  operator,  i,  is  also  used.  When  the 
inversion  operator  is  performed  on  a  molecule,  it’s  corresponding  wavefunction  is 
either  changed  or  unchanged.  Spectroscopists  keep  track  of  this  symmetry  property 
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by  the  use  of  the  subscripts  g  (gerade  or  even)  and  u  (ungerade  or  odd).  In  Dirac 
notation  [2], 


i\9)  =  Is)  (2.14) 

i\u)  =  — |  u)  (2-15) 

Therefore,  in  the  case  of  singlet  delta  oxygen  (a1A9)  the  subscript  g  refers  to  the 
fact  that  when  the  inversion  operator  is  applied,  it’s  original  wavefunction  remains 
unchanged. 

2-4-5  The  Herzberg  Diagram  and  Selection  Rules.  Quantum  mechanically, 
the  allowed  transitions  result  from  the  computation  of  the  following  equation. 

Qba  =  (^|Q|^a)  (2.16) 

Where  Qba  is  the  scalar  electric-quadrupole  moment,  Ta  and  dh  are  the  three-dimensional 
wavefunctions  for  the  a  and  b  states  respectively,  and  Q  is  the  electric-quadrupole  mo¬ 
ment  operator.  Of  course,  the  equation  doesn’t  cause  the  selection  rules,  but  more 
precisely  is  the  mathematical  model  of  what  occurs  physically  in  nature  between 
these  two  electronic  states.  To  avoid  performing  the  mathematics  every  time  one 
wanted  the  selection  rules,  Gerhard  Herzberg  developed  a  diagram  which  accounts 
for  these  allowed  transitions  (and  the  resulting  selection  rules)  and  symmetry  prop¬ 
erties  of  the  rotational  levels  for  a  given  transition  [8].  The  Herzberg  diagram  is  a 
pictorial  representation  that  neatly  encapsulates  all  of  the  selection  rules  and  dis¬ 
plays  the  allowed  transitions  between  rotational  levels  and  consequently  the  resulting 
lines  and  branches  that  are  detectable.  As  shown  in  Figure  2.10  five  branches  are 
allowed  for  the  b1^  <—  a1  Ag  system  of  02.  The  0,P,Q,R  and  S  branches  correspond 
to  AK  =  —2,  —1,0, 1,2  respectively4.  This  is  the  first  selection  rule  for  an  electric- 
quadrupole  transition.  Another  selection  rule  involves  the  electric-quadrupole  allowed 

4Recall  that  since  J  =  K,  then  AJ  =  AK  =  —2,  —1,0, 1,  2. 
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Figure  2.10:  This  is  the  Herzberg  diagram  for  the  <—  a1  Ag  transition.  This 
diagram  is  a  handy  tool  for  displaying  the  allowed  rotational  transitions  which  come 
about  from  the  various  selection  rules  for  this  particular  transition.  [8] 

transitions  between  electronic  states  of  02 •  This  selection  rule  is  AA  =  0,  ±1,  ±2.  For 
the  <—  cFAg  transition,  the  £  symbol  for  02(&)  means  A  =  0  and  the  A  symbol 
in  02(a)  means  A  =  2,  therefore,  AA  =  2.  The  other  selection  rules  for  an  electric 
quadrupole  transition  deal  with  the  molecular  symmetry  labels.  These  selection  rules 
are  g  g,  u  u,  +  +,  and  —  .  The  forbidden  transitions  are  g  **+  u  and 

+  «*•  -  [2]  [8]  [25], 

Referring  to  Figure  2.10,  it  can  be  seen  that  there  does  not  exist  any  transitions 
among  the  anti-symmetric  states.  The  reason  for  this  is  due  to  the  fact  that  the 
nuclear  spin  (I)  for  an  oxygen  atom  is  zero.  The  zero  quantity  then  classifies  the 
oxygen  atom,  and  consequently  the  molecule,  as  a  boson.  In  quantum  mechanics, 
the  total  wavefunction  for  a  bosonic  system  must  be  symmetric  according  to  Pauli’s 
Exclusion  Principle.  Looking  briefly  at  the  total  wavefunction  (ipTot),  we  have 

IpTot  'Pel  V'vib  Prof  Pnuc  (2.17) 

r 

Where  pei  is  the  electronic  wavefunction,  and  is  symmetric  for  b 1£^  <—  a1A9  due  to 
the  g  subscripts  for  each  state.  ipvn,  is  the  vibrational  wavefunction,  yet  this  is  solely 
dependent  on  the  internuclear  separation,  r,  so  this  is  also  symmetric.  As  mentioned 
earlier  I  =  0  for  oxygen,  so  the  nuclear  wavefunction  {ipnuc)  is  symmetric.  This  leaves 
the  rotational  wavefunction,  tprot .  However,  since  the  exclusion  principle  dictates  that 
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the  total  wavefunction  must  be  symmetric,  then  ip rot  must  also  be  symmetric.  This 
is  why  only  symmetric  transitions  occur  in  the  b x£+  <—  a1  Ag  transition. 

The  symmetry  properties  (i.e.  the  +  and  —  associated  with  each  electronic 
state)  depend  on  the  symmetry  properties  of  the  field  in  which  the  electrons  move. 
First,  consider  the  precession  of  the  L  vector  in  Figure  2.9.  If  an  angle  p  is  defined  in 
one  direction,  say  counterclockwise,  then  a  —  ip  will  exist  in  the  clockwise  direction. 
Therefore,  two  electronic  eigenfunctions  will  exist  [8], 


Ai  =  \e+i (2.18) 
ipei  =  xe~iAipi  (2.19) 

Since  the  Schrodinger  equation  is  linear,  the  linear  combination  of  the  above  two 

eigenfunctions  will  also  be  a  solution  to  the  same  eigenvalue.  The  total  electronic 
eigenfunction  can  be  combined  two  ways  [8], 

</,+  =  aXe+iAtpi  +  bXe~iAlfil  (2.20) 

^ei  =  aXe+lA<pl  ~  bXe~lAipl  (2.21) 

Where  a  and  b  are  arbitrary  constants.  Since  £  states  are  non-degenerate  (because 
A  =  0),  i\)ei  has  only  one  solution  and  is  no  longer  dependent  on  p\.  Therefore,  the 
eigenfunction  can  either  be  ip^  or  ip~v  In  fact  all  £  states  oscillate  between  having  a 
ip^  and  a  ip^,  however  the  labels  (+  or  — )  on  the  £  states  describe  which  electronic 
eigenfunction  state  they  begin  as.  For  example,  b x£+  will  start  out  being  symmetric 
and  consequently  all  positive  K  values  will  be  symmetric.  All  odd  K  values  are 
then  antisymmetric.  This  is  why  in  molecular  spectroscopy  all  £  states  are  labeled 
with  either  +  or  — .  The  state  a1Ag  however,  will  contain  both  the  symmetric  and 
the  antisymmetric  component  of  ipei  for  every  K  level.  This  is  why  II  states  and 
higher  have  a  two-fold  degeneracy.  All  of  this  is  shown  in  the  Herzberg  diagram  in 
Figure  2.10. 
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It  is  important  to  note  that  the  rotational  levels  of  a  diatomic  molecule  are 
classified  according  to  the  behavior  of  the  total  eigenfunction  (not  of  the  rotational 
eigenfunction  alone)  with  respect  to  reflection  at  the  origin.  A  rotational  level  is 
called  positive  or  negative  depending  on  whether  the  total  eigenfunction  remains 
unchanged  or  changes  sign  for  such  a  reflection.  This  property  is  also  called  parity 
and  corresponds  to  the  evenness  and  oddness  of  atomic  energy  levels.  [8] 

2.5  Spectral  Lineshape  Theory 

There  are  some  basic  concepts  of  spectroscopy  that  need  to  be  mentioned  for  the 
purpose  of  defining  commonly  used  phrases,  such  as  the  FWHM,  the  Lorentzian,  the 
Gaussian,  and  the  Voigt.  These  concepts  and  others  are  mentioned  in  the  following 
sections. 

2.5.1  Line  Widths.  Spectral  lines  are  never  completely  monochromatic  in 
absorption  or  emission  spectra.  If  they  were,  they  would  appear  as  infinitely  thin 
lines  at  their  center  frequency.  However,  there  is  always  a  spectral  distribution  I{y) 
around  the  central  frequency  r0  =  E\Ek  corresponding  to  a  molecular  transition 
with  an  energy  difference  A E=Ei~Ek  between  the  upper  and  the  lower  states.  The 
function  I{y )  in  the  vicinity  of  uq  is  called  the  line  profile,  and  the  frequency  interval 
Au=  \u2  —  v\\  is  called  the  Full  Width  at  Half  Maximum  (FWHM)  of  the  line  [4], 
There  is  another  line  width  that  is  used  called  the  Half  Width  at  Half  Maximum 
(HWHM).  The  HWHM  is  simply  half  the  FWHM.  Two  other  parts  of  the  line  profile 
are  the  line  kernel  and  the  line  wings.  The  kernel  is  simply  the  region  within  the 
FWHM  and  the  wings  are  the  region  outside  the  FWHM.  This  is  graphically  shown 
in  Figure  2.11  [4], 

If  the  isolated  molecule  or  atom  were  to  be  ideally  at  rest,  then  natural  lifetime 
broadening  would  be  the  dominate  feature  of  the  line  profile.  The  shape  of  this 
feature  is  called  the  Lorentzian  profile.  Classical  mechanics  is  used  to  describe  this 
frequency  distribution  as  a  damped  oscillator  and  ultimately  describes  the  shape  of 
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Figure  2.11:  This  is  a  general  line  profile.  [4] 
the  Lorentzian  profile5.  The  normalized  Lorentzian  line  function  is  defined  as  [19], 

Av 


(2.22) 


Where  An  is  the  FWHM5 6.  Pressure  broadening  also  contributes  to  the  Lorentzian 


lineshape  of  a  given  transition,  more... 

However,  real  molecules  are  never  at  rest.  There  will  always  be  some  thermal 
characteristic  to  the  molecule.  At  thermal  equilibrium,  the  molecules  of  a  gas  follow 
a  Maxwellian  velocity  distribution.  This  distribution,  which  is  fundamentally  Gaus¬ 
sian,  is  the  dominant  profile  for  the  lineshape.  Since  we  are  dealing  with  velocities, 
the  resulting  FWHM  is  often  called  the  Doppler  width.  This  width  is  the  major  con¬ 
tributor  to  the  lineshape  and  consequently  the  Lorentzian  lineshape  is  concealed  by 

5For  a  more  complete  explanation  refer  to  Laser  Spectroscopy  by  Derntroder.  [4] 

6Bernath  [2]  describes  the  FWHM  as  A  v  =  1/(27t rsp)  where  rsp  is  the  spontaneous  lifetime  of 
the  state.  Bernath  further  shows  that  the  FWHM  agrees  with  the  uncertainty  principle  A E  At  >  h , 


h 


1 

2  7 r 


(2.23) 


(2.24) 


20 


it.  Equation  2.25  is  the  form  of  the  Gaussian  line  function  [2]. 


(2.25) 


Where  m  is  the  mass  of  the  molecule,  c  is  the  speed  of  light  in  a  vacuum,  k  is 
Boltzmann’s  constant,  T  is  temperature  in  kelvin,  and  the  Doppler  FWHM  is  defined 
as,  [2]. 


(2.26) 


One  can  clearly  see  the  thermal  properties  of  the  molecule  contained  within  equa¬ 
tions  2.25  and  2.26. 

The  true  shape  of  the  spectral  line  is  neither  Lorentzian  nor  Gaussian,  but  a 
convolution  of  the  two,  called  a  Voigt  lineshape.  Equation  2.27  shows  the  result  of 


the  mathematics  behind  convoluting  equations  2.22  and  2.25  together  [19]. 


(2.27) 


Where  «o  is  the  area  under  the  Voigt  profile,  A vl  is  the  HWHM  of  the  Lorentzian 
and  Az/£>  is  the  Doppler  HWHM.  This  pressure  broadening  is  a  direct  consequence  of 
collisions  which  occur  due  to  the  molecules’  individual  velocities.  In  other  words,  the 
gas  is  pressure  broadening  itself.  Figure  2.12  compares  the  Voigt  lineshape  to  both 
the  Gaussian  and  the  Lorentzian  lineshapes. 

The  Voigt  profile  plays  a  crucial  role  in  the  analysis  of  this  experiment.  The 
raw  data  which  is  fit  against  a  Voigt  lineshape  is  used  to  determine  the  area  of  the 
line.  This  fit  area  is  necessary  in  calculating  the  peak  cross  section.  Also,  accurate 
measurements  of  the  Voigt  line  wings  enables  the  Lorentzian  FWHM  and  the  Gaussian 
FWHM  to  be  extracted.  Although  extracting  the  Lorentzian  and  Gaussian  lineshapes 
is  done  throughout  the  entire  experiment,  this  technique  is  especially  important  in 
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Figure  2.12:  The  Voigt  profile.  [4] 


determining  the  pressure  broadening  coefficients.  In  this  part  of  the  experiment  the 
Lorentzian  FWHM’s  of  a  line  are  plotted  against  the  pressure  they  were  scanned  at. 
The  slope  from  this  linear  relationship  is  precisely  the  pressure  broadening  coefficient 
of  that  particular  line.  It  is  interesting  to  note  that  at  low  pressures  (A ud  A/yp) 
the  Voigt  reduces  to  the  Gaussian  lineshape,  and  at  high  pressures  (A vd  -C  Aup)  the 
Voigt  reduces  to  the  Lorentzian  lineshape. 


2.5.2  Absorption  Parameters.  The  absorption  coefficient  relates  the  change 
in  flux  ( dF )  of  a  stream  of  photons  as  they  propagate  through  some  distance  (dx) 
containing  an  absorbing  species  to  their  original  flux  (F—  ^),  where  h  is  Planck’s 
constant  [2]  [9].  This  defines  the  absorption  coefficient  as, 


1  dF 
F  dx 


=  —a 


(2.28) 


The  rate  equation  for  a  two  level  system  must  be  considered  in  order  to  un¬ 
derstand  how  the  absorption  coefficient  ( a )  relates  to  the  integrated  absorption  cross 
section  (oy).  Since  this  experiment  deals  only  with  the  absorption  band, 

^  =  Bb^apNa  (2.29) 

Where  Bh^a  is  the  Einstein  coefficient  for  absorption,  p  is  the  radiation  density  (p  = 
AL),  and  both  Na  and  Nb  are  the  number  densities  of  02(a)  and  02(b)  respectively. 
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Substituting  the  definition  for  the  Einstein  B  coefficient  (namely  Bb^a  =  3^2  Qlo.gi.1'- 
Vba)) and  rearranging  terms,  it  can  be  shown  that  equation  2.29  becomes, 


dNb 

dt 


2?r2gg> 

3  e0hc 


g{v 


Vba)NaF 


(2.30) 


Where  Qba  is  the  electric-quadrupole  moment  that  governs  the  transition  between 
02(a)  and  02(b)1 .  The  symbol  v  is  the  frequency  of  the  observed  transition  and 
g(v  —  vba)  is  the  line  shape  function* 8.  The  integrated  absorption  cross  section  can 
now  be  defined  from  equation  2.30  as, 

97 t202 

ai(u)  =  -^ug(u-uba)  (2.31) 

o€qiIC 


This  is  the  ’’effective  area”  that  the  molecule  presents  to  that  stream  of  photons  of 
flux  F  [2], 

The  rate  at  which  02(a)  populates  02(b)  (i.e.^f)  is  proportional  to  the  change 
in  flux  of  the  photons  by  02(a)  over  some  distance,  and  can  be  expressed  as, 


dF  _  dNb 
dx  dt 


(2.32) 


Combining  equations  2.30,  2.31  and  2.32,  it  can  be  shown  that, 

dF 

—  =  -ajNbF  (2.33) 


Finally,  the  absorption  coefficient  can  be  shown  to  be  a  function  of  the  integrated 
absorption  cross  section  by  inserting  equation  2.33  into  equation  2.28, 


a  =  ajNb 


(2.34) 


'This  is  the  scalar  quantity  that  results  from  equation  2.16. 

8Which  is  typically  normalized  with  the  following  relation,  g(v)dv  =  1- 
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This  suggests  that  the  absorption  coefficient  of  an  absorbing  species  (in  this  case 
02(a))  is  ciirectly  proportional  to  the  effective  absorbing  area  of  that  absorbing  species 
multiplied  by  the  number  density  of  the  excited  state  it  is  populating  (in  this  case 
OoM). 

There  is  another  absorption  cross  section  (a)  that  is  frequently  reported  in 
papers,  and  this  is  the  cross  section  which  is  dependent  on  the  temperature  and 
pressure  of  the  optical  cavity.  The  units  for  this  cross  section  is  cm2.  This  value  is 
specific  to  a  certain  frequency.  So,  when  the  phrase  peak  cross  section  is  used,  this 
is  referring  to  the  cross  section  at  line  center,  vq.  This  cross  section  is  often  confused 
with  the  intrinsic  value,  called  the  integrated  absorption  cross  section  (  defined  in  this 

study  as  oy).  This  value  is  not  dependent  on  the  temperature  nor  the  pressure  and  is 

2 

reported  in  units  of  cm  , . 

2.5.3  Pressure  Broadening.  Pressure  broadening  refers  to  the  widening  of 
the  spectral  line  due  to  an  increase  in  pressure.  This  rise  in  pressure  increases  the 
probability  of  collisions  among  the  molecules  within  a  gaseous  mixture.  This  is  why 
pressure  broadening  is  often  referred  to  as  collisional  broadening.  To  conserve  energy 
and  momentum,  the  spectral  line’s  absorbency  (as  well  as  it’s  peak  cross  section) 
decreases  as  it  broadens.  This  shortening  and  widening  of  the  spectral  line  is  clearly 
shown  in  the  results  presented  in  Chapter  IV. 

In  previous  work  [21]  [22]  [19],  there  is  a  clear  dependence  on  the  rotational 
quantum  number  for  the  magnetic-dipole  allowed  transition  of  the  b <—  V3E~ 
system  of  oxygen.  Work  by  Steinfeld  et  al  [10]  reflects  the  views  of  many  experimental 
spectroscopists  that  rotationally-inelastic  collisions  are  the  dominant  mechanism  in 
molecular  pressure  broadening.  Therefore,  this  study  will  plot  this  trend  for  the  first 
time  for  the  Noxon  system. 
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III.  Experimental  Approach 


After  grasping  a  firm  understanding  of  the  theories  presented  in  Chapter  II,  it 
is  now  necessary  to  explain  the  experimental  approach  to  fully  appreciate  the 
data  obtained.  Therefore,  this  chapter  will  be  divided  into  three  main  sections.  The 
first  section  deals  with  the  technique  used  in  this  experiment.  The  second  section 
describes  the  apparatus,  while  the  third  section  explains  the  overall  procedure,  from 
optimizing  the  apparatus  to  analyzing  the  data. 

3.1  Multi-Pass  Absorption  Spectroscopy 

The  spectroscopic  method  used  in  this  experiment  is  aptly  named  multi-pass 
absorption  spectroscopy.  However,  the  name  is  really  a  general  phrase  for  the  methods 
employed  in  this  experiment.  It  is  called  multi-pass  because  the  cavity  has  the  ability 
to  bounce  light  on  the  order  of  80,000  times.  This  allows  a  gain  of  the  absorbance 
feature  to  occur  within  the  cavity.  The  word  absorption  refers  to  the  fact  that  this 
technique  deals  with  the  observation  of  absorption  lines  in  the  spectrum  as  opposed 
to  transmission  lines.  More  specifically,  however,  multi-pass  absorption  spectroscopy 
incorporates  cavity  ringdown  measurements  and  the  off-axis  integrated-cavity-output- 
spectroscopy  (ICOS)  technique.  These  methods  will  be  individually  discussed  below. 

3.1.1  The  Cavity  RingDown.  The  only  element  of  cavity  ringdown  spec¬ 
troscopy  (CRDS)  used  in  this  experiment  will  be  the  calculation  of  the  ringdown 
time.  The  ringdown  time,  r,  is  used  to  verify  how  well  the  laser  light  is  ’’trapped”  in 
the  optical  cavity.  In  order  to  acquire  the  best  performance  out  of  the  mirrors  in  the 
cavity,  they  must  be  devoid  of  any  additional  lossy  materials  such  as  water  and  dust 
particulates.  The  ringdown  time  will  therefore  give  a  measure  of  these  losses.  The 
value  of  the  ringdown  time  is  also  used  in  calculating  the  intensity  of  the  observed 
absorption  line.  This  will  be  further  explained  later  in  this  chapter. 

In  CRDS,  the  lifetime  of  a  photon  (This  is  simply  the  rate  at  which  it  takes 
light  to  enter  and  exit  the  cavity,  commonly  referred  to  as  the  ringdown  time)  in  an 
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optical  cavity  is  used  to  obtain  the  total  internal  cavity  losses  (per  pass)  as  a  function 
of  the  photon  frequency  [16].  When  these  losses  are  dominated  by  the  scatter  and 
transmission  of  the  two  mirrors  (i.e.  a  very  clean  cavity),  then  the  mirror  reflectivity 
can  be  obtained.  When  a  narrow  band  absorbing  species  is  present  (i.e.  02(a)), 
absolute  molecular  absorption  intensities  can  be  inferred  by  subtracting  the  baseline 
losses  of  the  cavity1  [16].  The  subtracting  of  this  baseline  will  be  further  explained  in 
the  Section  3.3. 


Figure  3.1:  The  1/e  point  on  the  downward  curve  of  this  square  wave  is  where  the 
ringdown  time  is  measured.  Although  the  x-axis  is  currently  in  point-space,  this  axis 
can  be  better  thought  of  as  a  time  axis  for  this  example. 

The  ringdown  time  is  acquired  in  this  experiment  by  first  modulating  the  laser 
signal  with  a  square  wave.  This  modulation  will  simulate  a  pulsed  laser  which  is 
necessary  in  measuring  the  ringdown  time  (As  shown  in  Figure  3.9).  At  the  end  of 
one  pulse,  the  output  waveform  has  a  rampdown  exponential  curve  which  corresponds 
to  the  time  decay  of  the  photon  population  leaving  the  cavity.  The  ringdown  time  is 
obtained  by  exponentially  fitting  this  time  decay  curve.  Typical  ringdown  times  are 
on  the  order  of  ca.  50 /as,  however,  this  experiment  has  acquired  ringdowns  as  high  as 
220 /as.  Optimizing  this  ringdown  time  is  crucial  when  detecting  weak  signals,  such 
as  those  generated  by  an  electric-quadrupole  transition.  To  understand  the  physical 

lrThis  statement  refers  to  the  relation  1  —  1-  See  a  further  explanation  in  Appendix  A. 
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nature  and  usefulness  of  the  ringdown  time,  it  is  necessary  to  develop  some  simple 
equations  which  will  relate  r  to  other  physical  quantities  of  the  cavity. 

Consider  when  the  reflectivity,  R ,  raised  to  the  number  of  trips  in  the  cavity, 
N,  is  equal  to  the  1/e  point  of  the  output  intensity2 


Rn  = 


1 

e 


Solve  for  N, 


N 


1 

h VR 


(3.1) 


(3.2) 


In  the  limit  of  high  reflectivity,  the  following  series  expansion  is  used, 

ln(x )  &  (x  —  1)  —  ^-(x  —  l)2  +  -(x  —  l)3  —  ...  (3.3) 

Zj  O 

Letting  R  =  x.  Therefore,  the  first  order  approximation  to  equation  (3.3)  becomes 


InR  «  R  —  1 


(3.4) 


Substituting  equation  (3.4)  into  equation  (3.2)  yields  the  following  relation  for  N, 


N 


1 


l-R 


(3.5) 


Now  a  relation  is  needed  for  the  ringdown  time,  r,  which  describes  how  long  it  takes 
light  to  travel  through  the  cavity  N  times,  at  a  total  distance  L  (the  length  of  the 


2The  ringdown  time  also  physically  refers  to  the  time  it  takes  the  output  intensity,  /,  to  equal 
the  incident  laser  intensity  divided  by  e  [14].  Therefore  the  ringdown  time  is  acquired  when  I  =  (" . 
This  relation  can  then  be  used  in  the  Beer-Lambert  formula  for  two  trips  (equation  (A. 3))  inside  the 
cavity,  I  =  IqR2  e~^aL .  Transmittancy  has  no  part  in  these  equations  because  we  are  only  concerned 
about  the  time  inside  the  cavity.  So,  —  =  IqR2 e~3aL .  Now,  Iq  cancels  out,  a  =  0  since  a  ringdown 
is  the  time  it  takes  light  to  travel  through  an  empty  cavity.  The  equation  then  becomes,  ^  =  R2 .  If 
we  then  consider  the  light  to  travel  through  the  cavity  N  times,  we  have  the  final  form  )  =  RN . 
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cavity).  Assuming  a  vacuum  state,  r  is  expressed  as, 


N  L 

T  =  - 

c 


where  c  is  the  speed  of  light  in  a  vacuum.  Solving  for  N, 


(3.6) 


N  = 


T  C 

~L 


(3.7) 


Finally,  set  equations  (3.5)  and  (3.7)  equal  and  solve  for  r, 


L 

c(l  -  R) 


(3.8) 


Equation  (3.8)  explicitly  shows  the  importance  of  high  reflectivity.  As  R  — »  1,  r  — >  oo. 
This  equation  will  help  the  experimenter  determine  if  the  laser,  mirrors,  detector,  etc. 
are  all  operating  properly,  since  maximum  r  is  desired.  Additionally,  this  equation 
will  allow  the  conversion  efficiency,  or  yield  of  02(a)  generated  by  the  plasma  to  be 
determined  for  a  given  data  set.  This  quantity  is  also  very  important  when  optimizing 
the  apparatus  for  maximum  02(a)  detection. 

Most  spectroscopists,  when  performing  CRDS  experiments,  use  a  high-powered 
pulsed  IR  laser.  A  pulsed  laser  is  the  obvious  choice  due  to  the  fact  that  measuring 
a  decay  rate  of  the  photon  population  is  at  the  heart  of  CRDS.  This  is  necessary, 
since  the  detector  needs  to  ’’see”  a  quanta  of  light  to  perform  the  required  decay  rate 
calculation  mentioned  above.  However,  pulsed  IR  sources  usually  do  not  have  the 
high  frequency  resolution  (narrow-band)  necessary  for  many  IR  applications,  such  as 
detecting  02(a).  Therefore,  tunable  diode  lasers  are  needed  to  provide  the  required 
narrow-band  resolution.  However,  diode  lasers  cannot  simply  be  used  in  place  of  the 
pulsed  laser  in  a  CRDS  set-up.  One  reason  is  the  fact  that  CRDS  set-ups  require  laser 
light  to  be  injected  into  the  cavity  on  the  optical  axis.  Although  this  provides  a  stable 
cavity  configuration,  the  noise  level  associated  with  the  optical  feedback  to  the  diode 
laser  due  to  the  on-axis  configuration  is  drastic  and  counter-productive.  Therefore, 


an  off-axis  design  is  required.  Specifically,  the  design  developed  by  O’Keefe  [16]  will 
be  utilized  in  this  experiment.  This  design  is  capable  of  tuning  over  many  modes. 

3.1.2  Off-Axis  ICOS.  The  theory  of  off-axis  paths  in  optical  cavities  have 
been  around  since  the  1960 ’s.  The  reflections  of  light  rays  in  such  configurations  are 
spatially  separated  until  the  ’’re-entrant”  condition  is  fulfilled  [12].  The  ’’re-entrant” 
condition  refers  to  when  the  light  ray  begins  to  retrace  it’s  path  within  the  cavity. 
This  spatial  separation  greatly  reduces  the  optical  feedback  for  the  tunable  laser  diode. 
Depending  on  the  design  of  the  cavity,  the  re-entrant  condition  can  be  met  from  a 
few  passes  to  an  infinite  number  of  passes.  In  order  to  achieve  cavity  enhancement  of 
the  absorption  feature,  our  experiment  requires  a  large  number  of  passes  before  the 
re-entrant  condition  is  met.  This  is  accomplished  with  the  use  of  parabolic  mirrors 
and  made  use  of  their  natural  astigmatism.  Specifically  the  mirrors  are  designed  after 
the  parabolic  mirrors  in  the  Herriott  cell.  When  astigmatic  mirrors  are  used  the  laser 
light  creates  a  Lissajous  spot  pattern  on  the  mirrors  and  throughout  the  cavity.  This 
pattern  alleviates  the  central  spot  pattern  along  the  optical  axis  which  creates  the 
unwanted  optical  feedback  to  the  tunable  laser  diode.  Because  of  this,  off-axis  ICOS  is 
truly  the  only  method  currently  developed  to  utilize  the  diode  lasers  with  the  required 
resolution  to  detect  02(a). 

Once  the  re-entrant  condition  is  met  with  the  desired  amount  of  passes,  the 
next  factor  to  be  concerned  about  is  the  ratio  between  the  cavity  free  spectral  range 
(FSR)  and  the  bandwidth  of  the  absorbing  feature.  In  order  to  scan  linearly  through 
the  cavity  with  the  laser  frequency,  the  FSR  of  the  cavity  must  be  narrower  than 
the  absorption  feature.  (Due  to  the  incredibly  long  path  length  of  our  design,  the 
modes  within  the  cavity  basically  overlap  and  we  have  effectively  a  continuous  mode. 
This  allows  for  the  laser  to  scan  linearly  very  effectively.  This  linear  relationship  is 
crucial  in  the  analysis  of  the  data,  since  a  polynomial  is  fit  to  the  baseline  and  is 
subtracted  from  the  absorption  feature.  The  more  the  baseline  resembles  a  line,  the 
more  accurate  the  analysis. 
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Starting  from  Beers  Law  and  incorporating  the  above  concepts  we  arrive  at  the 


following  equation  which  is  the  heart  of  off-axis  ICOS  theory, 


A  /  _  GA 

I(aL  =  0)  ““  1  +  GA 


(3.9) 


Where  G=  AC  is  the  cavity  enhancement  factor,  A—  1  —  e~aL  is  the  single-pass 
absorption,  and  A I—I(aL  =  0 )—I(aL)  is  the  change  in  the  cavity  output.  Essen¬ 
tially,  this  equation  says  that  the  absorption  information  (A)  is  contained  within  that 
change  in  the  cavity  output  [16] .  This  allows  for  a  direct  measurement  of  the  absorp¬ 
tion  feature.  A  more  detailed  explanation  of  each  mathematical  step  is  presented  in 
Appendix  A. 


3. 2  The  Apparatus 

The  apparatus  used  in  this  experiment  is  designed  for  the  off-axis  integrated- 
cavity-output-spectroscopy  (ICOS)  technique,  and  utilizes  the  method  of  cavity  ring- 
down  to  determine  the  decay  time,  (r),  of  the  light  within  the  cavity.  The  reason 
why  this  particular  device  is  unique  is  because  of  it’s  ability  to  ’’trap”  light  to  such 
a  high  degree.  In  Section  2.1  it  was  explained  that  off-axis  ICOS  was  developed  by 
Anthony  O’Keefe  in  2001.  However,  O’Keefe  was  only  able  to  achieve  approximately 
400  passes  in  his  cavity,  whereas  this  device  allows  for  greater  than  80,000  passes. 

Since  the  resolving  power  of  the  apparatus  depends  on  the  number  of  passes 
the  photons  will  complete  in  the  cavity,  it  becomes  critical  to  have  highly  reflective 
mirrors.  The  mirrors  used  in  this  experiment  are  coated  to  highly  reflect  the  near  IR 
(1400  —  1600  nm)  and  have  a  radius  of  curvature  of  ca.  lm.  Reflectivities  of  more  than 
99.9987%  were  achieved  and  confirmed  with  the  ringdown  time3.  As  mentioned  above, 
the  shape  of  the  mirrors  were  modeled  after  the  Herriott  cell  (parabolic  mirrors)  which 
is  designed  to  ’’trap”  light.  Therefore,  the  heart  of  this  experiment  lies  within  the  two 

3Use  equation  (3.8)  and  solve  for  R.  In  this  experiment  the  highest  ringdown  time  achieved  was 
220 /j,s.  This  corresponds  to  R  =  99.9988% 
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mirrors  placed  in  the  cavity.  Keeping  the  mirrors  devoid  of  any  particulates,  especially 
water,  becomes  an  ongoing  concern  whenever  a  series  of  data  is  taken.  After  cleaning 
the  mirrors,  they  stayed  clean  for  about  10  days. 

To  generate  the  needed  02(a)  species  of  oxygen,  a  mixture  of  ultra  high  purity 
helium  (99.999%)  and  research  grade  oxygen  is  used.  This  mixture  is  then  sent  to 
a  quartz  tube  via  mass  flow  meters  which  control  both  the  flow  rates  and  the  ratio 
between  the  two  gases.  A  McCarroll  cavity  is  placed  over  a  region  of  the  quartz 
tube  so  that  a  microwave  discharge  can  be  applied  across  the  mixture.  A  Tesla 
Coil  is  used  to  apply  the  required  initial  ionization  to  ignite  the  gas  to  create  the 
plasma.  The  plasma  will  then  sustain  itself  and  no  further  voltage  is  required.  The 
McCarroll  cavity  is  cooled  with  pressurized  air  to  prevent  over-heating.  This  plasma 
will  generate  the  singlet  delta  oxygen  needed  for  the  experiment.  The  02(a)  then 
travels  through  the  main  cavity  of  the  apparatus  and  is  evacuated  using  a  vacuum 
pump.  A  schematic  diagram  of  the  experiment  setup  is  shown  below  in  Figure  3.2. 
The  main  cavity  contains  the  highly  reflective  parabolic  mirrors  at  either  end  to 


Figure  3.2:  This  is  the  experiment  setup.  The  length  of  the  optical  cavity  is  82 

cm,  the  microwave  power  supply  operated  at  80  W,  and  the  average  pressure  in  the 
cavity  was  ca.  23  torr. 
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increase  the  path  length  of  the  transmitted  laser  light.  A  tunable  laser  diode,  which 
is  coupled  to  a  fiber,  is  first  connected  to  a  light  attenuator,  which  splits  the  light  into 
two  separate  fibers.  This  first  fiber  has  a  1:0.99  light  intensity  ratio,  and  is  mounted 
near  the  entrance  of  the  cavity.  The  light  beam  enters  the  cavity  through  the  back 
of  one  of  the  mirrors.  The  light  then  exits  the  cavity  by  means  of  the  extremely  low 
transmittancy  of  the  back  mirror,  passes  through  a  lens  which  focuses  the  photons 
onto  the  InGaAs  detector  which  then  sends  the  data  to  a  computer  where  it  is  stored 
as  l’s  and  0’s  in  a  text  hie.  The  second  fiber  on  the  attenuator  has  a  10:1  light 
intensity  ratio.  This  decrease  in  intensity  is  needed  for  the  wavemeter  which  is  used 
to  determine  the  absolute  location  of  the  spectral  lines.  The  wavemeter  specified  a 
safe  range  of  light  intensities  from  0  to  20  dBm.  After  going  through  the  attenuator, 
the  laser  beam  decreased  from  90  dBm  to  about  10  dBm. 

3.2.1  The  Diode  Lasers.  Five  distributed  feedback  (DFB)  laser  diodes  were 
utilized  in  this  experiment.  Tunable  DFB  laser  diodes  were  preferred  for  their  reso¬ 
lution  with  the  added  convenience  of  being  able  to  change  their  scanning  wavelength 
peak  by  controlling  the  amount  of  current  being  sent  to  it.  Broad  tuning  is  achieved 
by  adjusting  the  temperature  of  the  diode,  and  fine  tuning  is  achieved  by  adjust¬ 
ing  the  current  supplied  to  the  diode.  These  lasers  were  designed  to  scan  at  and 
around  the  following  wavelengths:  1494.48  nm,  1499.71  nm,  1505.74 nm,  1512.575  nm, 
and  1515.63  nm.  Diode  lasers  are  excellent  for  this  type  of  work  since  their  scanning 
wavelengths  are  controlled  via  a  thermal  coupler.  This  experiment  will  operate  in  a 
safe  temperature  range  of  10°C  to  49°C.  With  this  operating  range,  the  lasers  can 
only  search  over  a  wavelength  range  of  4  nm.  Furthermore,  each  laser  has  a  linewidth 
of  0.0001  crrT1.  This  aids  in  the  resolution  required  to  observe  ultra-weak  transitions 
which  had  linewidths  of  ca.  0.008  cm-1. 

An  etalon  measurement  of  the  diode  laser  was  made  during  each  scan  of  a  line. 
This  is  necessary  in  forming  a  baseline  during  the  analysis,  since  the  etalon  allows  for 
the  conversion  of  each  data  point  into  frequency-space. 
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3. 3  Procedure 


As  with  any  experiment,  the  procedure  for  this  experiment  began  with  calibra¬ 
tion  and  cleaning.  In  the  first  section  below,  a  brief  description  is  given  as  to  how 
the  tunable  laser  diodes  were  calibrated.  This  calibration  was  crucial  in  determining 
the  exact  locations  of  the  absolute  line  positions.  The  next  step  in  the  procedure  was 
to  optimize  various  elements  of  the  equipment  to  produce  the  best  results.  For  this 
experiment,  the  best  results  correspond  to  high  signal  strength  (i.e.  maximize  the 
ringdown  time)  and  low  noise.  Therefore,  a  brief  history  is  laid  out  which  describes 
the  optimization  process  for  this  experiment.  Next,  a  section  is  provided  below  that 
briefly  discusses  how  the  data  was  taken.  Finally,  the  last  section  provides  a  short 
description  of  how  the  data  was  analyzed.  This  is  necessary,  because  the  results  pre¬ 
sented  in  both  Chapter  IV  and  Appendix  C  are  actually  the  end  result  of  a  series 
of  steps  performed  by  a  computer  program  called  IGOR  Pro®.  IGOR  Pro®  is  a 
powerful  data  analysis  program. 

3.3.1  Calibration  of  Lasers.  The  purpose  of  calibrating  the  DFB  laser  diode 
is  to  extract  the  required  temperature  setting  which  corresponds  to  a  line  in  the  02(a) 
spectra.  Since  these  laser  diodes  only  scan  through  0.33  cm'1,  only  one  line  can  be 
seen  at  a  time.  This  narrow  scanning  window  is  another  reason  why  calibrating  the 
lasers  is  so  important. 

This  experiment  utilized  a  wavemeter  to  calibrate  the  lasers.  As  mentioned 
above,  a  safe  operating  range  of  10 °C  to  49°C  was  used  in  this  experiment.  Anything 
below  10 °C  risked  condensation  forming  on  the  diode  and  anything  above  49°C  risked 
the  chance  of  thermally  damaging  the  electronics  within  the  diode.  Starting  at  10°C, 
the  laser  was  fed  into  the  wavemeter  which  read  the  peak  wavelength  of  the  laser  in 
wavemeters  (cm-1).  The  laser  temperature  was  then  increased  in  steps  of  3°C'  until 
25°C',  and  then  in  steps  of  2 °C  until  49°C'.  At  each  step,  the  peak  wavelength  was 
recorded.  A  graph  was  then  generated  for  each  laser  which  showed  it’s  temperature  as 
a  function  of  frequency  (in  cm-1).  See  Figure  3.3  for  an  example  of  one  such  graph.  As 
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The  1499.71  nm  Laser  Diode 
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Figure  3.3:  Temperature  dependence  for  the  1499.71  nm  laser 

Table  3.1:  Lines  that  are  within  the  safe  temperature  range  for  the 
1499.71  nm  tunable  laser  diode. 


Line 

Energy  (cm 

Laser 

Temperature  (°C) 

S2 

6662.169255 

29.5888 

R3 

6653.662739 

46.5181 

R5 

6658.343271 

37.2031 

R7 

6662.661218 

28.6097 

R9 

6666.615169 

20.7407 

Rll 

6670.203603 

13.5992 

expected,  there  is  a  clear  linear  relationship  between  the  temperature  of  the  diode  and 
the  peak  (or  center)  wavelength  it  will  scan  around.  A  linear  regression  is  performed 
on  the  data  in  order  to  extract  the  slope  and  y-intercept.  These  parameters  were  then 
used  to  calculate  the  temperature  the  laser  needed  to  be  at  for  a  given  spectral  line’s 
energy.  For  example,  for  the  laser  shown  in  Figure  3.3,  Table  3.2  lists  the  lines  that 
can  be  observed  with  their  corresponding  laser  temperatures.  Tables  such  as  these 
were  generated  for  each  of  the  five  laser  diodes. 

3.3.2  Optimization. 
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Table  3.2:  List  of  the  slopes  and  y-intercepts  obtained  for  the  five 

lasers. 


Laser  (ran) 

Slope  ( °C/cm  4) 

Y-intercept  (°C) 

1494.71 

1499.71 
1505.74 
1512.575 

1515.63 

-2.0905  ±  0.0046 

-1.9902  ±  0.0003 

-2.0956  ±  0.0053 

-2.1323  ±  0.0037 
-2.1345  ±  0.0038 

14008  ±  31.0 

13288  ±  1.89 

13942  ±  35.5 

14114  ±  24.2 

14109  ±  24.9 

3.3.2. 1  Flow  Rate  Dependence.  One  of  the  most  important  elements  is 
understanding  the  relationship  between  the  concentration  of  02(a)  and  the  flow  rate 
ratio  of  oxygen  and  helium.  I11  order  to  accomplish  this  we  first  set  the  laser  to  the 
strongest  line  in  the  Noxon  Band,  Q(8).  Fixing  the  helium  flow  at  4000  seem4,  the 
oxygen  was  first  brought  in  at  40  seem,  then  100  seem,  200  seem,  400  seem,  and  finally 
500  seem.  At  each  increment  of  oxygen  the  Q(8)  line  was  recorded  and  then  analyzed 
using  IGOR  Pro  ®.  From  the  analysis,  the  number  density  (Na)  of  02(a)5  ,  the  yield 


4sccm  stands  for  Standard  Cubic  Centimeter  per  Minute.  Where  ’’Standard”  refers  to  Standard 
Temperature  (273  K)  and  Pressure  (1  atm)  or  STP. 

5The  number  density  is  found  by  the  following  method.  Since  the  single-pass  absorption  is 
defined  as  Ax  =  1  —  e~axNxL (where  o\  is  the  cross  section  of  a  certain  wavelength,  of  02(a),  Nx 
is  the  number  density  of  a  certain  wavelength  of  02(a),  and  L  is  the  length  of  the  cavity),  then  we 
can  say  axNxL  =  ln(l  —  ^4a)-  From  our  analysis  we  measure  both  G  (from  the  ringdown  time)  and 
GA  (which  is  the  peak  signal  strength).  Therefore,  we  are  able  to  calculate  Ax  from  the  relation 
^4.  Now  we  can  write  <j\ N\L  =  ln(l  —  ^=4).  Next  we  need  two  relations  for  the  number  density  of 
02(a),  namely  [02(a)]  =  NTfb  (where  fb  is  the  boltzmann  factor  for  a  certain  wavelength  and  Nt 
is  the  total  number  density  of  02(a)),  and  A\  =  \02(a)]Lax.  Setting  each  equal  to  [02(a)]  and  then 
each  other,  we  generate  a  formula  for  the  total  number  density  of  02(a), 


Nt 


Ax 

C\fbL 


(3.10) 


If  we  series  expand  the  exponent  in  A\  =  l~e~axNxL,  we  find  that  the  first  order  approximation  will 
be  Ax  =  <jxNxL.  This  means  that  Ax  ss  ln(  1  —  qA).  Inputting  this  result  for  Ax  into  equation  3.10, 
the  total  number  density  of  02(a)  can  finally  be  written  as, 


ln(  1  —  ^4) 

Nt  =  — - g-L 

<yxfbL 


(3.11) 
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of  02(a) 6,  and  the  relative  intensity  (GA)  were  all  able  to  be  extracted.  Figures  3.4 
and  3.5  illustrate  the  relation  among  these  quantities.  Referring  to  Figure  3.4,  we  see 


Figure  3.4:  Flow  rate  dependence  for  singlet  delta  oxygen 


Figure  3.5:  Relationship  between  the  intensity  of  the  signal  and  the  number  density 
of  02(a). 


that  the  number  density  of  02(a)  increases  to  a  point  as  more  oxygen  is  introduced. 
However,  as  more  oxygen  is  introduced  the  yield  of  02(a)  decreases.  This  is  mainly 

6The  yield  is  obtained  by  simply  dividing  the  number  density  of  02  (a)  by  the  total  number  density 
of  oxygen. 
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due  to  the  self-quenching  mechanism  of  02(o). 


02(a)  +  02(d)  — >  02(X)  +  02(X)  (3.12) 

Therefore,  with  more  02(a)  produced  there  is  a  higher  probability  of  self-quenching 
among  the  02(a).  Although  having  a  background  in  kinetic  theory  would  have  pro¬ 
vided  an  a  priori  to  this  part  of  the  experiment,  there  was  not  enough  time  to  explore 
this  subject  in  detail.  Trial  and  error  proved  to  fair  just  as  well,  however,  and  it 
was  concluded  that  setting  the  oxygen  flow  meter  to  400  seem  (10%  of  the  helium) 
provided  the  higher  signal  strength  (Figure  3.5),  even  though  the  yield  was  not  the 
highest  it  could  have  been.  For  measuring  the  exact  location  of  the  lines  in  02(a), 
it  is  more  important  to  maximize  the  signal-to- noise  ratio  than  it  is  to  maximize  the 
yield. 


3. 3. 2. 2  Power  Dependence.  Another  aspect  to  optimize  in  this  exper¬ 
iment  is  the  amount  of  power  the  microwave  power  supply  generates.  As  mentioned 
in  Section  3.2,  microwaves  are  sent  across  a  quartz  tube  via  a  McCarroll  cavity.  A 
Tesla  coil  is  used  to  ignite  the  plasma,  but  afterwards  it  becomes  self  sustaining  and 
stable  for  several  hours.  The  microwave  power  supply  has  two  analog  readouts,  one 
for  the  forward  power  and  the  other  for  the  reflected  power.  By  adjusting  controls  on 
the  McCarroll  cavity,  it  was  possible  to  reduce  the  reflected  power  to  about  zero. 

Once  again  the  Q(8)  line  was  chosen  as  the  test  line.  The  line  was  first  analyzed 
at  60  Watts  forward,  and  then  stepped  in  increments  of  10  Watts  until  90  Watts. 
From  the  data  collected  and  from  prior  experience  (Williams  et  al  previous  work),  it 
was  concluded  to  use  80  Watts  forward  power. 

3. 3. 2. 3  Alignment.  There  are  two  very  critical  areas  that  require 
alignment.  The  first  is  the  alignment  of  the  cavity.  Since  this  apparatus  is  designed 
for  off-axis  ICOS  it  is  extremely  important  that  the  mirrors  were  aligned  to  give  the 
maximum  ringdown  time.  Recall  from  above  that  the  ringdown  time  refers  to  the 
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amount  of  time  it  takes  a  photon  to  leave  the  cavity.  Also  recall  that  as  the  ringdown 
time  increases,  the  cavity  enhancement  factor  increases  as  well.  So,  a  great  deal  of 
time  was  devoted  to  the  proper  alignment  of  this  cavity.  A  computer  was  used  to 
show  the  ringdown  time  in  real-time,  so  that  the  ringdown  time  could  be  seen  with 
each  adjustment  made.  This  alignment  was  only  done  at  the  start  of  the  experiment. 

The  second  item  that  required  careful  alignment  was  the  laser  itself.  The  tunable 
laser  diodes  were  coupled  to  a  fiber  at  one  end  and  the  other  end  was  mounted  near  the 
entrance  of  the  cavity,  ft  is  this  mount  that  controls  the  direction  of  the  laser  beam. 
Once  again  the  ringdown  time  was  monitored  as  the  mount  was  minutely  adjusted. 
This  particular  alignment  required  a  bit  of  give  and  take  between  maximizing  the 
ringdown  time  and  minimizing  the  noise.  The  ringdown  time  would  increase  as  the 
laser  beam  approached  the  optical  axis  (but  not  on  the  axis),  but  the  noise  level 
would  increase  because  the  optical  feedback  would  increase  as  the  beam  approached 
the  optical  axis.  This  alignment  was  probably  checked  every  other  week  for  about 
three  months. 

3. 3. 2. 4  Various  Other  Optimizations.  There  are  other  elements  that 
required  optimizing  on  a  more  regular  basis.  One  is  the  lens  at  the  end  of  the  cavity. 
Prior  to  collecting  a  series  of  lines  for  a  given  laser,  the  lens  was  adjusted  to  make 
sure  the  light  exiting  the  cavity  was  focused  directly  onto  the  InGaAs  detector.  This 
ensured  maximum  intensity.  Another  aspect  of  optimizing  is  cleaning  the  mirrors.  In 
earlier  sections  it  was  emphasized  that  the  reflectivity  of  the  mirrors  play  a  paramount 
role  in  the  cavity  enhancement  factor.  This  being  said,  the  mirrors  were  cleaned  from 
once  a  week  to  a  few  times  a  day,  depending  on  the  value  of  the  ringdown  time.  The 
mirror  reflectivities  would  degrade  mainly  due  to  water  build  up.  This  experiment 
used  only  helium  and  oxygen,  but  there  was  still  traces  of  water  from  the  ambient 
air.  Even  though  a  vacuum  pump  was  used  to  pump  out  the  ambient  air  from  the 
cavity,  an  ideal  vacuum  can  never  be  obtained.  As  a  a  result,  there  will  always  be 
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some  traces  of  water  in  the  cavity7.  In  fact,  water  lines  were  frequently  observed  when 
the  diode  lasers  were  dialed  through  a  range  of  energies.  There  were  even  cases  when 
a  water  line  either  overlapped  with  the  02(a)  line  (such  as  R(5),  P(ll)  and  S(6))  or 
was  within  the  0.33  cm window  near  the  02(a)  line. 

3. 3. 2. 5  Optimizing  for  Pressure  Broadening.  After  the  lines  were 
located,  the  next  step  to  the  experiment  is  to  pressure  broaden  select  lines  in  the  02(a) 
spectra.  Now  the  importance  is  more  on  the  yield  of  02(a)  than  on  the  intensity.  This 
is  due  to  the  fact  that  at  higher  pressures  the  probability  of  self-quenching  increases 
and  it  becomes  difficult  to  keep  02(a)  sustained  in  the  cavity.  For  example,  if  the  flow 
rate  of  oxygen  was  kept  at  400  seem  and  the  total  pressure  in  the  cavity  is  increased 
to  100  Torr  (from  about  20  Torr),  then  the  line  would  simply  vanish  into  the  noise 
level.  I11  fact,  under  these  conditions  the  Q(8)  line  was  difficult  to  see  at  50  Torr. 
Having  this  knowledge,  and  the  information  from  Figure  3.4,  the  oxygen  was  chosen 
to  be  set  at  50  seem.  This  corresponds  to  about  1%  of  the  helium. 

In  order  to  extract  a  pressure  broadening  coefficient  for  a  given  line,  the  Lorentzian 
linewidth  must  be  determined  from  the  fit  data.  It  becomes  critical  then  that  the 
shape  of  the  Voigt  lineshape,  especially  the  wings,  is  not  distorted  as  the  laser  scans 
through  the  line.  For  determining  the  positions  of  the  lines  a  scanning  rate  of  25  Hz 
was  more  than  adequate.  For  pressure  broadening,  however,  this  was  too  fast.  Many 
trial  runs  were  conducted  using  25  Hz,  10  Hz,  and  5  Hz  scan  rates.  The  distortion 
can  be  observed  when  the  Voigt  profile  is  fit  to  the  data.  Specifically,  the  raw  data 
would  lean  to  the  left  or  right  on  either  side  of  the  fitted  Voigt  profile8.  The  25 
Hz  scan  proved  too  fast,  because  the  Voigt  fit  clearly  misses  the  data  in  Figure  3.6. 
However,  there  was  virtually  no  distortion  observed  in  either  the  10  or  the  5  Hz  scan 
(Figures  3.7  and  3.8  respectively).  In  fact,  there  didn’t  appear  to  be  any  significant 
difference  between  the  two  other  than  the  5  Hz  scan  produced  more  noise.  Actually, 

7The  porous  nature  of  the  plastic  gas  lines  used  in  the  apparatus  allowed  water  to  always  be 
present. 

sFor  examples  of  this  refer  to  the  data  presented  in  Appendix  C. 
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Frequency  (cm ') 

Figure  3.6:  This  is  the  Q8  line  scanned  at  25  Hz.  Total  pressure  in  the  cavity  was 
ca.  25  torr,  with  a  ringdown  time  of  203  [a s.  The  data  was  averaged  10000  times  and 
binned  40  times. 


Frequency  (cm ') 

Figure  3.7:  This  is  the  Q8  line  scanned  at  10  Hz.  Total  pressure  in  the  cavity  was 
ca.  25  torr,  with  a  ringdown  time  of  203  /as.  The  data  was  averaged  4000  times  and 
binned  100  times. 
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Frequency  (cm  ) 


Figure  3.8:  This  is  the  Q8  line  scanned  at  5  Hz.  Total  pressure  in  the  cavity  was 
ca.  25  torr,  with  a  ringdown  time  of  203  /is.  The  data  was  averaged  1000  times  and 
binned  200  times. 


the  reason  why  the  5  Hz  scan  appeared  to  produce  more  noise  was  because  the  25 
Hz  scan  was  averaged  10,000  times,  the  10  Hz  scan  4,000  times  and  the  5  Hz  scan 
1,000  times.  These  averages  were  chosen  in  order  to  keep  each  scan  per  line  under  ten 
minutes.  When  pressure  broadening,  each  line  was  scanned  at  four  different  pressures 
form  25  Torr  to  100  Torr,  so  each  line  took  about  45  minutes  to  scan. 


3.3.3  Taking  Data.  Section  3.2  dealt  with  the  apparatus  in  the  sense  of 
how  the  singlet  delta  oxygen  is  produced  and  how  it  is  collected  by  the  laser  and  the 
detector.  This  section  will  instead  focus  on  the  computer  program  used  to  collect  the 
raw  data. 

There  are  actually  two  programs  used  when  taking  the  data.  The  first  is  called 
Noxon.  This  program  displays  and  exponentially  fits  the  ringdown  curve  in  real-time 
and  calculates  the  ringdown  time.  This  program  was  essential  in  both  determining 
the  ringdown  time  for  a  given  series  of  scans,  and  aided  in  the  maintenance  of  the 
apparatus  since  a  low  ringdown  time  corresponded  to  dirty  mirrors  or  bad  alignment. 
Figure  3.9  shows  an  example  of  the  ringdown  curve.  The  shape  is  a  result  from 
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Figure  3.9:  The  ringdown  curve  for  r  =  219  fis. 

modulating  the  laser  with  a  square  wave  which  simulates  the  pulsed  lasers  mentioned 
earlier.  When  the  laser  is  ’’switched  off’  the  exponential  decay  of  the  light  within 
the  cavity  is  clearly  seen.  The  light  will  also  decay  ”up”  until  the  cavity  reaches 
steady-state  conditions,  which  corresponds  to  when  the  re-entrant  condition  is  met. 
This  occurs  at  approximately  80,000  passes  in  the  cavity  for  this  experiment. 

The  other  computer  program  is  aptly  called  I  COS.  This  program  is  the  main 
interface  between  the  user  and  the  apparatus.  This  program  displays  both  the  real¬ 
time  signal  and  the  averaged  signal  on  the  same  screen.  It  is  this  program  that  allows 
the  user  to  adjust  the  sampling  rate,  the  repetition  rate,  the  signal  type,  the  voltage 
offset  and  amplitude,  the  input  channel,  and  the  number  of  averages.  The  repetition 
rate  (or  rep  rate)  means  that  each  data  point  will  be  sampled  a  certain  number  of 
times;  this  is  measured  in  Hertz.  Whereas  the  sampling  rate  is  measured  in  points  per 
second  and  refers  to  the  total  amount  of  samples  taken  by  the  detector,  regardless  of 
repetitions.  For  example  if  the  sampling  rate  is  500,000  Hz  and  the  rep  rate  is  25  Hz, 
then  there  will  be  20,000  data  points  that  will  describe  a  line.  When  a  ringdown  is 
measured,  using  this  program,  the  rep  rate  is  set  at  100  Hz,  so  there  will  be  5,000  data 
points  to  describe  the  ringdown  curve.  This  is  shown  along  the  x-axis  in  Figure  ??. 
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The  signal  type  refers  to  the  type  of  waveform  with  which  the  laser  can  be  modulated. 
A  square  waveform  is  used  for  generating  a  ringdown  curve  and  a  sawtooth  waveform 
is  used  when  scanning  a  line.  The  input  channel  allows  the  user  to  toggle  between 
the  InGaAs  detector  and  the  etalon  measurement.  Finally  the  number  of  averages 
refers  to  how  many  times  a  frequency  will  be  averaged.  When  determining  the  line’s 
position  the  number  of  averages  was  set  to  10,000  times.  It  was  mentioned  above  that 
at  25  Hz  rep  rate  there  will  be  20,000  data  points  to  describe  the  line,  so  each  of  the 
20,000  points  will  now  be  averaged  10,000  times.  This  is  why  it  took  approximately 
ten  minutes  to  scan  each  line.  The  following  section  will  cover  the  process  of  analyzing 
the  data  that  these  programs  saved  as  text  hies. 

3.3.4  Analyzing  Data.  The  raw  data  obtained  from  the  InGaAs  detector, 
which  was  saved  as  a  text  hie,  is  then  processed  using  another  computer  program 
called  IGOR  Pro  ®.  Williams  wrote  macro  program  hies  specifically  for  the  analysis 
of  this  type  of  data. 

When  the  data  is  hrst  uploaded  into  the  program  it  is  convenient  to  overlay 
the  raw  data  trace  (the  absorption  line  in  the  sawtooth  waveform)  with  the  etalon 
function.  This  is  displayed  in  Figure  3.10. 

Explain  better  (include  the  binning  process) 

The  etalon  function  tells  the  user  how  many  fringes  were  present  within  the 
laser’s  0.33  cm-1  scanning  ’’window”  for  each  line.  More  importantly  however,  the 
etalon  provides  the  function  which  will  determine  the  non-linear  relationship  between 
point-space  (labeled  as  Point  Number)  and  frequency-space.  Although  it  is  not  ap¬ 
parent  in  this  example  (Figure  3.10),  the  fringe  peaks  are  not  equally  spaced  in  point 
space.  In  fact,  if  the  laser  was  allowed  to  scan  at  a  higher  rep  rate,  then  many 
fringes  would  be  seen  and  there  would  be  a  visible  bunching  of  fringe  peaks  of  varying 
separation  intervals.  However,  an  etalon  will  alleviate  this  problem  since  the  free- 
spectral  range9  of  that  etalon  is  the  spacing  between  the  fringes  in  frequency-space. 

9The  free-spectral  range  is  2  GHz  for  the  etalon  used. 
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Figure  3.10:  The  raw  data  for  R(7)  with  the  etalon  function. 


Frequency  (cm ') 
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Figure  3.11:  The  etalon  is  used  to  convert  the  R(7)  data  to  frequency-space. 
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Therefore,  the  step  being  performed  between  Figure  3.10  and  Figure  3.11  is  first  iden¬ 
tifying  the  mathematical  function  which  will  transform  the  data  points  to  points  in 
frequency-space,  and  then  performing  this  function  on  the  raw  data.  The  result  of  this 
transformation  is  shown  in  Figure  3.11  and  the  new  absorption  feature  is  now  plotted 
as  a  function  of  some  arbitrary  frequency  of  units  cm-1.  The  data  in  frequency-space 
is  bowed  out  from  the  original  data  because  the  non-linear  transformation  function  is 
parabolic. 

After  the  data  has  been  corrected  for  frequency-space,  the  data  is  then  binned. 
Binning  refers  to  the  process  of  averaging  a  group  of  points  together  to  form  one  point. 
This  averaging  will  also  reduce  the  statistical  noise  level.  It  was  mentioned  earlier 
that  the  sampling  rate  and  scan  rate  both  play  a  factor  as  to  how  many  points  will 
represent  a  particular  line  in  the  spectra.  Since  the  sampling  rate  was  fixed  at  500 
kHz  throughout  the  entire  experiment,  a  10  Hz  scan  produced  50,000  points  and  a  25 
Hz  scan  produced  20,000  points.  Now,  in  this  experiment  10  Hz  scans  were  binned 
at  100  times,  which  reduced  50,000  points  to  500  points,  and  the  25  Hz  scans  were 
binned  at  40  times,  which  reduced  20,000  points  to  500  points.  These  two  scan  rates 
(10  Hz  and  25  Hz)  can  now  be  properly  compared  to  one  another,  given  the  fact  that 
they  now  have  the  same  number  of  data  points  representing  the  absorption  feature. 

The  next  step  is  where  physics  plays  a  role.  ICOS  theory  provides  a  relationship 
among  the  intensity  of  the  empty  cavity  (. I(aL  =  0)),  the  intensity  of  the  cavity  with 
the  absorber  present  ( I{aL ))  and  the  Voigt  profile  of  that  absorber  {GA).  These 
quantities  are  related  to  each  through  the  following  equation10. 


I(aL  =  0) 
I{ptL) 


GA 


(3.13) 


This  is  the  exact  form  used  in  the  macro  hie  in  IGOR  Pro  ®.  How  this  relation  is 
used  is  better  understood  if  some  direct  comparisons  are  made  to  (3.13).  Since  G 

lnA  full  mathematical  development  of  this  equation  from  the  Beer-Lambert  formula  is  provided 
in  Appendix  A. 
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is  the  cavity  enhancement  factor  and  A  is  the  single-pass  absorption  (SPA)  feature, 
then  GA  can  be  viewed  as  the  enhanced  single-pass  absorption  of  the  absorption 
species;  in  other  words  G(SPA).  Since  I(aL  =  0)  is  the  output  intensity  of  the 
empty  cavity,  then  we  can  equate  this  quantity  to  the  baseline  polynomial  fit  (/ (aL  = 
0)  =  Baseline).  Since  I(aL)  is  the  output  intensity  of  the  cavity  with  the  absorbing 
species  present,  then  we  can  equate  this  quantity  to  the  raw  data  taken  from  the 
InGaAs  detector  (. I(aL )  =  Rawdata ).  Writing  this  out  in  equation  form  we  have, 


Baseline 

RawData 


1  =  G(SPA) 


(3.14) 


Figure  3.12  shows  the  polynomial  fit  to  the  baseline,  the  Voigt  fit  to  the  peak  and 


Frequency  (cm ') 

0.05  0.00  -0.05  -0.10  -0.15  -0.20 


Figure  3.12:  The  absorbance  is  calculated. 


the  corresponding  residuals.  According  to  equation  (3.14)  the  macro  then  divides  the 
baseline  out  of  the  raw  data,  subtracts  one  (this  will  invert  the  absorption  peak)  and 
then  calculates  the  value  of  GA.  This  GA  value  will  correspond  to  the  peak  of  the 
line  shape.  The  result  is  shown  in  Figure  3.13.  Finally,  the  last  macro  is  performed 
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and  those  result  are  shown  in  Figure  3.14.  Three  main  calculations  take  place  in  this 
macro.  First,  another  more  accurate  Voigt  profile  is  fit  to  the  data.  This  is  done 
now  that  the  baseline  is  removed,  so  the  values  for  the  Gaussian  and  the  Lorentzian 
lineshape  functions  will  be  more  accurate.  The  next  calculation  involves  the  cross 
section  of  the  line.  The  third  and  final  calculation  finds  the  peak  location  in  frequency- 
space  by  inserting  the  value  obtained  from  the  wavemeter.  The  wavemeter  was  used 
during  the  collection  of  each  line,  via  the  fiber  optic  attenuator.  The  wavemeter 
provided  the  frequency  of  the  peak  of  the  laser  line  width  to  an  accuracy  of  0.004 
cm,-1.  The  peak  of  the  laser  line  width  correspond  to  the  center  point  of  the  raw  data 
in  point-space  (i.e.  the  center  of  the  graph  in  Figure  3.10).  Since  the  line  is  already 
transformed  to  the  laser’s  frequency-space,  the  program  simply  replaces  the  center 
frequency  with  the  peak  frequency  found  by  the  wavemeter.  The  x-axis  is  then  scaled 
to  this  new  frequency,  and  the  result  is  shown  in  Figure  3.14. 

3.3.4- 1  Taking  Data  for  Pressure  Broadening.  When  acquiring  data 
for  pressure  broadening  all  of  the  above  procedures  are  still  followed.  The  only  differ- 
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Figure  3.14:  A  Voigt  profile  is  fit  to  the  peak  and  the  peak  cross  section  is  calculated. 
This  is  the  final  step  in  the  analysis  process. 


ence  is  the  fact  that  each  line  must  now  be  scanned  four  times  in  a  row,  corresponding 
to  each  of  the  four  pressures.  Those  four  pressures  are  25  Torr,  50  Torr,  75  Torr  and 
100  Torr.  These  four  pressures  are  unique  in  and  of  themselves  since  most  spectro- 
scopists  start  at  or  above  100  Torr. 
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IV.  Results  and  Discussion 


The  results  presented  in  this  chapter  prove  with  great  success  the  ability  off- 
axis  ICOS  has  to  detect  extremely  weak  absorption  signals.  All  results  are 
focused  on  the  (1,0)  band  of  the  b <—  a1 2  Ag  system  in  oxygen.  The  first  set  of 
results  confirms  the  absolute  line  positions  of  21  lines  and  are  compared  with  the 
expected  positions  predicted  from  the  Cheah  and  Lee  results  [27]  of  the  02(a)  and 
02(b)  term  energies.  Thirteen  of  these  lines  have  never  been  observed  before.  These 
lines  encompass  all  five  branches,  which  are  predicted  by  quantum  mechanics  and  are 
shown  in  the  Herzberg  diagram  in  Section  2.2.5.  The  next  set  of  results  verify  the 
integrated  cross  sections1  using  a  form  of  Boltzmann  analysis.  The  final  set  of  results 
list  pressure  broadening  coefficients  for  select  lines  in  the  (1,0)  band  and  verifies, 
for  the  first  time,  the  rotational  dependence  on  pressure  broadening  for  an  electric- 
quadrupole  transition.  Although  the  results  are  mainly  presented  in  tabular  form  in 
this  chapter,  all  of  the  analyzed  lines  are  presented  graphically  in  Appendix  C. 

4-1  The  Absolute  Line  Positions 

The  direct  quantitative  determination  of  the  absolute  line  positions  for  the 
0,P,R  and  S  branches  of  the  b <—  a1  Ag  (1,0)  transition  of  oxygen  had  not  been 
accomplished  before  this  experiment.  Williams  et  al  [24],  however,  already  measured 
the  Q  branch  quantitatively  from  J=2  through  J=16.  The  same  Q’s  were  measured 
in  this  experiment  and  the  values  agree  to  within  0.001cm-1  of  Williams’  values. 

All  of  the  measured  values  from  this  experiment,  presented  in  Table  4.1,  are 
compared  to  their  predicted  values  calculated  from  the  Cheah  and  Lee  results2,3. 
Prior  to  the  invention  of  the  off-axis  ICOS  technique,  the  energies  for  the  Noxon  band 
could  only  be  calculated  from  subtracting  the  term  energies  a1  Ag  <—  X3£-  from 
b1^  <—  X3T,~.  The  measured  values  differ  from  the  predicted  values  by  0.0001cm-1 

1These  integrated  cross  sections  were  first  calculated  by  Dr.  Skip  Williams 

2The  process  of  how  Cheah  and  Lee’s  results  were  used  to  determine  the  band  energies  are 
presented  in  Appendix  C,  in  the  form  of  a  Matlab®  code. 

3The  spectroscopic  values  used  for  the  Hamiltonian  were  taken  from  the  Cheah  and  Lee  paper  [27]. 
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Table  4.1:  Absolute  Line  Positions  in  Wavenumbers  (cm  4). 


Line 

Quantitative 

Measurement 

lcr  Error 

Predicted  from 
Cheah  et  al 

Difference 

06 

6611.1418 

±0.0046 

6611.13567869935 

0.0061 

08 

6598.8226 

±0.0046 

6598.81569124817 

0.0069 

P9 

6614.4826 

±0.0046 

6614.47987116784 

0.0027 

P13 

6599.3977 

±0.0046 

6599.39161739536 

0.0060 

P15 

6591.3228 

±0.0046 

6591.31499534225 

0.0077 

Q2 

6642.9513 

±0.0047 

6642.94970151615 

0.0015 

Q4 

6642.3227 

±0.0046 

6642.32129012862 

0.0013 

Q6 

6641.3324 

±0.0046 

6641.33355116894 

0.0011 

Q8 

6639.9878 

±0.0046 

6639.98617083616 

0.0016 

Q10 

6638.2786 

±0.0046 

6638.27872121984 

0.0001 

Q12 

6636.2118 

±0.0046 

6636.21066030018 

0.0011 

Q14 

6633.7802 

±0.0046 

6633.78133194786 

0.0011 

Q16 

6630.9987 

±0.0046 

6630.98996592416 

0.0087 

R3 

6653.6667 

±0.0047 

6653.66273921216 

0.0039 

R7 

6662.6576 

±0.0047 

6662.6612178316 

0.0035 

R9 

6666.6174 

±0.0047 

6666.61516949024 

0.0022 

Rll 

6670.2048 

±0.0047 

6670.20360286338 

0.0012 

R13 

6673.4308 

±0.0047 

6673.42488042498 

0.0059 

R17 

6678.7650 

±0.0047 

6678.75884746335 

0.0061 

R19 

6680.8693 

±0.0047 

6680.86769134183 

0.0016 

S8 

6692.1316 

±0.0047 

6692.12146915857 

0.0101 

for  Q(10),  to  0.01cm”1  for  S(8).  The  S(8)  line  was  the  weakest  line  measured  {GA  = 
2.47554  x  10”3  for  S(8)  compared  to  GA  =  1.99853  x  10-2  for  Q(8))  in  this  experiment 
and  was  the  most  difficult  to  fit.  It  is  clearly  shown  in  Appendix  C  that  the  fitted 
peak  slightly  misses  the  ’real’  peak  of  the  data.  This  was  the  best  fit,  however,  for 

this  line. 

The  error  has  been  derived  from  the  wavelength  accuracy  of  the  wavemeter  used 
to  locate  the  peak  center  of  the  laser  linewidth  when  scanning  through  an  absorption 
line.  For  example,  the  calibration  chart  that  came  with  the  wavemeter  indicated 
that  for  a  line  at  1500nm  the  wavemeter  would  have  an  error  (A  A)  of  1.05  ppm4. 

4parts  per  million 
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Therefore,  the  error  in  wavenumbers  was  found  using  the  following  equation. 


A  A 

\v\  = 

A2 

The  measured  results  in  Table  4.1  agree  well  with  the  known  spectroscopic  theory. 
This  means  that  the  off-axis  ICOS  technique  does  indeed  find  extremely  weak  transi¬ 
tions,  such  as  produced  from  electric-quadrupole  transition. 

The  Q(8)  line  shown  in  Figure  4.1  is  an  example  of  what  the  line  looks  like  in 
frequency-space  after  being  fully  analyzed  from  the  raw  data  collected  by  the  appara¬ 
tus5.  The  line’s  cavity  enhanced  absorption  ( GA )  is  plotted  as  a  function  of  frequency 
(cm^1).  The  black  line  represents  the  actual  data  which  has  been  transformed  non¬ 
linear  ly  and  has  had  it’s  baseline  removed.  The  red  line  is  the  Voight  fit  to  the  data. 
Recall  that  the  peak  cross  section,  which  is  along  the  right  y-axis,  can  be  thought 
of  as  a  cavity  enhanced  integrated  absorption  cross  section.  The  residuals  above  the 
peak  represent  the  difference  between  the  raw  data  and  the  Voight  fit.  The  best 
signal  to  noise  ratio  (SNR)  obtained  in  this  experiment  was  40: l6.  Although  much 
time  and  effort  was  devoted  to  increasing  this  ratio,  it  still  fell  short  of  the  150:1 
results  Williams  et  al.  [24]  had  achieved.  Williams  was  able  to  achieve  such  a  high 
SNR  because  he  had  both  a  higher  ringdown  time  (250 /is  compared  to  220 /as)  and 
a  higher  yield  of  02(a)  (5%  conversion  efficiency  compared  to  2.7%).  However,  this 
experiment  did  achieve  noise  levels  lower  than  Williams  et  al. 

4- 2  Verification  of  the  Integrated  Absorption  Cross  Sections 

This  experiment  also  successfully  verified  the  intrinsic  integrated  absorption 
cross  sections  for  all  the  lines  of  Table  4.1.  This  was  done  by  using  a  form  of  Boltz¬ 
mann  analysis.  These  integrated  absorption  cross  sections  were  neither  measured  nor 
calculated  in  this  experiment,  but  were  first  calculated  by  Williams  et.al.  [24], 

5See  Chapter  III  for  a  detailed  explanation  of  the  procedure. 

6This  ratio  was  calculated  for  the  Q(8)  line.  The  worst  SNR  obtained  was  5:1  for  the  very  weak 
S(8)  line. 
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Frequency  (cm ') 

Figure  4.1:  The  Q(8)  line  position  at  6639.9878cm-1.  The  total  pressure  in  the 
cavity  was  ca.  25  torr,  the  ringdown  time  was  216  /is,  and  the  rep  rate  was  25  Hz. 
The  data  was  binned  50  times. 

4-2.1  Boltzmann  Analysis.  To  verify  the  integrated  absorption  cross  sections 
listed  in  Table  4.2,  the  following  three  equations  were  used. 

fB  =  (2S  +  1)(2J  +  l)e~%*  (4.2) 


[02(a)]  =  NafB  (4.3) 


A  =  [02{a)\LaI  (4.4) 


Where  fB  is  the  boltzmann  factor,  Ej  is  the  rotational  energy,  T  is  the  rotational 
temperature,  [02(a)]  is  the  concentration  of  02(a),  Na  is  the  number  density  of  02(a), 
A  is  the  peak  area  from  the  Voight  fit7,  L  is  the  length  of  the  cavity,  and  oy  is  the 
integrated  absorption  cross  section.  After  substitution  and  manipulation  of  the  above 

7This  is  found  from  the  following  relation,  ^4. 
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2 

Table  4.2:  Verified  Integrated  Cross  Sections  (  cm_ , ) . 


Line 

Integrated  Absorption 

Cross  Sections,  07 

06 

1.50233753287255x  10~25 

08 

1.34644197315733x  10-25 

P9 

4.44880627547526  xl0~25 

P13 

4.23990875270618xl0”25 

P15 

4.16985296122742xl0”25 

Q2 

4.44797314593381  xl0“25 

Q4 

5.45732707987711  xl0“25 

Q6 

5.65692593379x  10-25 

Q8 

5.72787157770784x  10“25 

Q10 

5.75900341017908x  10”25 

Q12 

5.77335389737773x  10~25 

Q14 

5.77909588510437  x  10”25 

Q16 

5.77973206313335x  10~25 

R3 

2.23476375652766  xl0~25 

R7 

3.14138016279576xl0-25 

R9 

3.31260656007808xl0-25 

Rll 

3.42615296599517  x  10~25 

R13 

3.5073662886619xl0-25 

R17 

3.61623411380298xl0-25 

R19 

3.65446514511864x  10~25 

S8 

6.8986524783927  x  10-26 

equations  it  can  be  easily  shown  that, 

ln  (o-/(25  +  l)(2J+l))  =  ~kf  ^ 

If  the  lines  are  plotted  with  their  In  ^  ai(2S+i)(2 J+i))  vafues  as  a  function  of  their 
individual  rotational  energy,  Ej,  the  slope  that  results  will  simply  be  —  yy-  From 
this  slope,  the  temperature  within  the  cavity  is  obtained.  Now,  when  the  O 2(a) 
leaves  the  plasma  and  travels  to  the  cavity  the  temperature  of  the  02(a)  will  decrease 
quickly  to  a  value  near  room  temperature  (~  300/1 ).  Therefore,  the  slope  must 
yield  a  temperature  near  room  temperature  in  order  to  verify  the  integrated  cross 
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sections.  All  21  lines  have  been  used  in  this  experiment  to  generate  the  slope  shown 


Figure  4.2:  Boltzmann  analysis. 

in  Figure  4.2.  A  temperature  of  314/1  ±  7 K  was  obtained  from  the  slope  of  the  linear 
regression  of  the  data.  This  result  compares  extremely  well  to  published  results. 
In  2004,  Williams  et  al.  [24]  calculated  a  temperature  of  312/i  ±  14 K  using  the 
same  Boltzmann  analysis  on  the  same  transition  band  of  oxygen8.  In  1985,  Fink 
et.al.  [6]  calculated  a  temperature  of  313/1  ±  10 K  using  the  same  Boltzmann  analysis 
from  the  emission  spectra  of  the  (0,0)  band  of  the  >  a1  Ag  system  of  oxygen. 

More  importantly,  however,  this  temperature  of  314/1  ±  7 K  is  consistent  with  the 
integrated  absorption  cross  sections  calculated  by  Williams.  Therefore,  the  calculated 
integrated  cross  sections  correlating  to  the  lines  in  this  experiment  have  been  verified 
experimentally. 

This  temperature  plays  a  crucial  role  in  extracting  the  Lorentzian  line  widths 
of  the  lines  when  pressure  broadened.  This  is  done  by  fixing  the  Gaussian  width  at 
this  temperature.  This  can  be  done  since  the  Gaussian  profile  is  a  result  of  Doppler 
broadening  which  is  dependent  on  the  temperature.  The  following  equation  illustrates 

86W+  <-  (1,0) 
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this  relationship. 


l2kTln(2 ) 

Aud  —  2uq\  - g —  (4.6) 

V  m 

Where  A  vo  is  the  Doppler  FWHM,  vq  is  the  natural  FWHM,  T  is  the  temperature, 
and  m  is  the  mass  of  the  molecule  or  atom.  The  pressure  broadening  results  will  be 
presented  in  the  following  section 


4-3  Pressure  Dependence  of  Line  Width  (FWHM) 


The  final  goal  for  this  experiment  is  to  pressure  broaden  select  lines  in  order 
to  extract  a  pressure  broadening  coefficient.  The  ten  strongest  lines  were  chosen 
for  their  ability  to  be  detectable  at  100  Torr9.  Four  pressures  were  used  in  this 
experiment;  25,  50,  75,  and  100  Torr.  Figures  4. 3-4. 6  graphically  illustrates  the 
Q(8)  line  being  pressure  broadened  at  four  pressures.  To  obtain  the  pressure 


1 .0  — , 
?  0.5 

°  0.0 
*  -0.5  - 
-1.0  J 


-1-1-L.l-l-l  I  I  I  I  I  l-l.Ll l.l-l  I  I  I  I  I  I  I  I  I  I  1 


V  "X 


broadening  coefficient,  the  resulting  Lorentzian  linewidth  are  then  plotted  against 


9The  cavity  was  typically  at  20  Torr  with  the  back  valve  completely  open 
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Figure  4.5:  The  Q(8)  line  pressure  broadened  at  75  Torr 
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^ WVV V V  m\f — W  WV 

pressure,  shown  in  Figure  4.7.  Each  data  point  is  a  full  ten  minute  scan  at  the 
specified  pressure.  The  error  bars  are  a  2cr  uncertainty  in  the  fit  to  the  Lorenztian 
FWHM.  To  acquire  decent  survey  data  of  the  Q(8)  line  sixteen  scans  were  taken, 
four  for  each  of  the  four  mentioned  pressures.  After  plotting  the  results,  a  linear 
regression  was  taken  with  a  weighted  error.  The  parameters  from  this  line  equation 
yielded  two  important  quantities.  The  most  important  parameter  was  the  slope,  since 
this  is  equal  to  the  peak’s  pressure  broadening  coefficient,  Once  the  the  pressure 
broadening  coefficient  is  determined,  the  Lorentzian  FWHM  (zq/2)  of  the  line  can 
then  be  calculated  at  any  given  pressure  (p)  from  the  relation  [2], 


^1/2  =  7 pP  +  intercept 


(4.7) 


The  pressure  broadening  coefficients  for  the  chosen  ten  lines  are  presented  in  Table  4.3 
along  with  their  ler  uncertainties. 

The  second  parameter  is  the  y- intercept.  This  quantity  proved  to  be  both 
interesting  and  perplexing.  In  the  absence  of  pressure,  equation  4.7  states  that  the 
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Pressure  (Torr) 

Figure  4.7:  Slope  of  regression  line  is  the  pressure  broadening  coefficient  in 

cm-1  Torr-1. 

FWHM  should  be  zero.  This  means  the  linear  regression  line  in  Figure  ??  needs  to 
pass  through  the  origin  oe  near  the  origin  within  reasonable  uncertainty.  The  data 
in  this  experiment  did  not  show  this.  One  cause  for  this  y-intercept  could  be  the 
signal  response  function  of  the  instrument  used  in  this  experiment.  This  is  possible 
because  this  experiment  resolves  linewidths  as  small  as  0.001  cm-1.  Further  work  will 
be  required  to  fully  explain  this  issue.  The  y-intercepts  along  with  their  error  is  also 
presented  in  Table  4.3. 

4-3.1  Rotational  (J)  Dependence  on  Pressure  Broadening  Coefficients.  The 
final  graph  of  this  chapter  (Figure  4.8)  depicts  the  rotational  dependence  of  the  pres¬ 
sure  broadening  coefficients.  Newman  and  Orr-Ewing  [26]  created  a  similar  graph  for 
the  a1  Ag  <—  (0,0)  transition.  The  trend  for  their  data  and  the  trend  shown 

in  Figure  4.8  both  suggest  that  as  the  rotational  quantum  number  increases,  the  in¬ 
fluence  that  the  collisions  have  on  the  Lorentzian  lineshape  decreases.  Pope,  for  his 
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Table  4.3:  Pressure  Broadening  Coefficients. 


Line 

Pressure 
Broadening 
Coefficient 
(cmT1  /torr) 

la  Error 
[cm-1  /torr) 

Y-lntercept 

(cm'1) 

lcr  Error 
(cm”1) 

P9 

4.83xlO”05 

±2.76xlO”06 

4.37xlO"04 

±1.54xlO”04 

Q4 

5.81xlO”05 

±3.51xlO”06 

6.42  xlO"04 

±2.30xl0”°4 

Q6 

5.06xl0”°5 

±2.72xlO”06 

4.75  xlO”04 

±1.55xlO”04 

Q8 

5.34xlO”05 

±2.43  xlO~06 

4.01  xlO”04 

±1.33xlO”04 

Q10 

4.98x  10”05 

±5.51  xlO-06 

3.18xlO”04 

±2.97x  10”04 

Q12 

4.19xlO”05 

±2.91  xlO-06 

7.99  xlO”04 

±1.65xlO”04 

Q14 

3.76x  10”05 

±6.59  xlO"06 

7.75  xlO”04 

±3.84xlO”04 

Q16 

3.61xlO”05 

±6.87xlO”06 

8.95xlO”04 

±3.66xlO”04 

R7 

5.77xlO”05 

±5.17xlO"06 

2.48  xlO”04 

±2.93xlO”04 

R9 

4.07xl0”°5 

±5.67xlO~06 

8.03xl0”°4 

±3.11xlO”04 

dissertation  [19]  [21]  [22],  plotted  the  same  trend  for  the  A  band  of  oxygen10.  This 
experiment  documents  for  the  first  time  a  rotational  dependence  on  the  Lorentzian 
lineshape  for  an  elcctric-quadrupolc  transition.  These  results  are  consistent  with  the 
known  theory  that  rotationally-inelastic  collisions  are  the  dominant  mechanism  in 
molecular  pressure  broadening  [19]  [21]  [22], 


10The  A  band  refers  to  the  A3£+  <—  X3£g  transition  in  oxygen.  This  transition  is  governed  by 
the  electric  dipole  moment. 
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Pressure  Broadening  Coefficient, y  (cm  /  torr) 
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Figure  4.8:  Rotational  dependence  on  the  pressure  broadening  coefficients. 


60 


V.  Summary  and  Recommendations 


All  of  the  main  researdi  goals  were  achieved  through  the  course  of  this  exper¬ 
iment.  The  results  presented  in  Chapter  IV  were  acquired  and  compared  to 
theory  with  great  success  and  accuracy.  In  addition  to  this,  the  results  also  establish 
off-axis  ICOS  as  a  proven  spectroscopic  technique  to  observe  ultra-weak  transitions, 
such  as  those  exclusively  dependent  on  the  electric-quadrupole  moment.  In  this  ex¬ 
periment,  twenty-one  lines  were  found,  their  integrated  absorption  cross  sections  were 
verified,  and  ten  lines  were  pressure  broadened  with  helium.  Time  did  not  allow, 
however,  the  ability  to  pressure  broaden  the  ten  lines  with  other  noble  gases,  such  as 
Ne,  Ar,  Kr  or  Xe.  This  chapter  will  provide  a  short  summary  of  the  results  and  then 
provide  recommendations  that  will  aid  in  the  continuation  of  this  or  related  research. 

5. 1  Summary 

The  results  of  Chapter  IV  verifies  the  sensitivity  of  the  spectroscopic  method 
used  in  this  experiment  to  locate  the  <—  a1  Ag  (1,0)  band  energies.  Absolute  line 

positions  were  determined  to  within  0.004cm_1.  The  measured  absolute  line  positions 
of  the  twenty-one  lines  differed  from  their  predicted  values  by  about  0.005  cm-1,  which 
is  on  the  order  of  the  la  error  from  the  wavemeter  (See  Table  4.1).  The  degree  of 
accuracy  achieved  in  this  research  effort  may  result  in  a  slight  correction  to  existing 
known  spectroscopic  constants  for  the  02(a)  and  02(b)  term  energies,  in  addition  to 
calculating  the  spectroscopic  constants  for  the  b 1X+  <—  a1  Ag  (1,0)  band  itself,  which 
has  not  been  done  before.  The  accuracy  of  these  line  positions  may  also  benefit  the 
HITRAN  database1. 

The  rotational  temperature  of  314  ±  7K ,  derived  from  a  Boltzmann  analysis  of 
the  21  detected  lines,  was  also  found  to  be  consistent  with  the  integrated  absorption 
cross  sections  calculated  by  Williams  et  al.  This  rotational  temperature  also  compares 
nicely  with  similarly  published  results.  Williams  et  al  calculated  312  ±  14 K  for  the 

1This  database  catalogs  all  known  absorption  lines  of  the  atmosphere  and  is  maintained  by  MIT. 
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same  band,  and  Fink  et  al  calculated  313  ±  10/i  for  the  emission  of  the  a1  A g 

(0,0)  band  [6]  [24], 

Pressure  broadening  coefficients  for  ten  lines  were  obtained  ranging  from  7 p  = 
3.61  x  10 /torr  to  5.81  x  10-6cm-1/forr.  These  coefficients  are  simply  the 
slope  of  a  regression  line  through  the  data  set  for  each  particular  line.  Since  these 
specific  lines  had  not  been  pressure  broadened  before,  there  were  no  published  val¬ 
ues  to  compare  these  results  to.  The  values  for  the  coefficients  were  reasonable  and 
decreased  as  the  rotational  quantum  number  increased  from  J=4  to  J=16.  This  par¬ 
ticular  property  of  pressure  broadening  has  been  observed  and  reported  for  other 
transitions  [26]  [19]  [21]  [22],  Being  able  to  graph  this  trend  was  the  only  method 
used  to  determine  if  the  coefficients  collected  in  this  experiment  were  indeed  reason¬ 
able.  Duplicating  the  pressure  broadening  analysis  done  in  this  research  effort  will  be 
necessary  in  solidifying  the  results  found. 

In  terms  of  applications,  this  research  effort  clearly  proved  the  ability  of  off- 
axis  ICOS  to  detect  singlet  delta  oxygen,  which  is  of  great  concern  to  the  COIL 
community.  From  it’s  inception  in  the  70’s  to  present-day  technology,  there  has  not 
been  a  method  yet  to  detect  the  concentration  of  02(a)  generated  from  the  SOG  of 
a  COIL.  This  apparatus  works  so  well  that  the  Chemical  Branch  at  Kirtland  AFB 
has  just  acquired  a  system  similar  to  the  apparatus  used  in  this  research  effort.  The 
spectroscopic  data  documented  in  this  thesis  should  be  a  great  foundation  for  their 
research  and  integration  efforts  with  the  COIL. 

The  atmospheric  community  will  most  benefit  from  the  pressure  broadening 
analysis.  The  coefficients  derived  in  this  study  should  aid  in  the  understanding  of  the 
radiation  transfer  of  the  upper  atmosphere.  Now  that  the  Lorentzian  widths  can  be 
determined  at  any  pressure  (in  torr),  these  coefficients  will  provide  a  better  model  for 
absorption  processes  which  take  place  in  the  upper  atmosphere. 
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5.2  Recommendations  For  Further  Research 


There  are  areas  of  this  experiment  that  have  the  ability  to  be  either  expanded 
upon  or  added  to  for  the  purpose  of  creating  a  new  thesis  research  topic.  Some  of 
those  potential  topics  are  presented  below. 

One  idea  is  to  acquire  an  RF  discharge  to  generate  the  plasma  instead  of  the 
microwave  discharge  that  was  used.  This  may  be  worth  exploring  for  the  purpose  of 
increasing  the  yield  of  02(a).  It  was  mentioned  in  previous  sections  that  Williams 
et.al.  [24]  obtained  conversion  efficiencies  of  5%,  which  was  double  of  what  was  gener¬ 
ated  in  this  experiment.  Part  of  the  reason  may  have  been  the  higher  ringdown  time, 
but  another  reason  could  be  the  RF  discharge  that  they  used.  With  higher  yields, 
some  of  the  more  weaker  lines  may  become  strong  enough  to  analyze  and  report  on. 

Another  topic  involves  the  area  of  pressure  broadening.  This  area  was  only 
slightly  explored  in  this  experiment  since  time  only  allowed  for  the  use  of  He.  In 
fact,  some  lines  pressure  broadened  in  this  experiment  by  He  were  only  scanned  four 
or  eight  times,  whereas  the  Q(8)  line  was  scanned  sixteen  times.  Therefore,  better 
survey  data  could  be  taken  of  all  the  lines  in  this  experiment  using  helium. 

Another  option  would  be  to  use  the  off-axis  ICOS  apparatus  to  analyze  a  com¬ 
pletely  different  ultra-weak  transition  besides  the  Noxon  Band.  This  will  require  the 
ordering  of  new  tunable  diode  lasers,  gases  and  mirrors  that  will  reflect  the  new  de¬ 
sired  frequency  range  of  interest.  If  new  mirrors  are  ordered,  then  funds  on  the  order 
of  $10,  000  must  be  considered  since  each  mirror  costs  approximately  $5,  000. 
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Appendix  A.  Mathematical  Development  of  GA 


A.l  Beers  Law  to  Off-Axis  ICOS  Theory 


According  to  Beers  Law,  the  output  intensity  with  an  absorbing  species  present 
( I(aL ))  is  related  to  the  incident  laser  intensity  (Jo)  through  the  following 
equation, 


I{aL)  =  I0e~a  L  (A.l) 

where  a  is  the  absorption  coefficient  of  the  molecular  absorber  (a  =  a(No  —  Ni)), 
and  L  is  the  length  of  the  main  cavity.  For  a  single  pass  in  the  main  cavity,  the 
transmittancy  of  the  mirrors  need  to  be  accounted  for. 

I(aL)  =  I0T2e~aL  (A.2) 

Where  T  is  the  transmission  coefficient  for  each  mirror  (in  our  case  the  mirrors  are 
identical).  However,  since  the  bMf  <—  a1Ag  system  is  such  a  weak  transition,  it  would 
be  advantageous  to  increase  the  output  intensity.  One  way  to  do  this  is  to  design  a 
multi-pass  off-axis  cavity.  The  benefit  is  more  apparent  through  the  development 
of  the  following  equations.  Equation  (A.2)  is  the  relation  for  light  simply  traveling 
through  the  cavity.  If  the  light  were  to  reflect  off  both  mirrors  (also  known  as  one 
roundtrip)  before  leaving  the  cavity,  then  it  would  have  traveled  three  times  the  length 
of  the  cavity.  This  is  expressed  in  the  following  relation, 

I(aL )  =  I0T2R2e~3aL  (A.3) 

where  R  is  the  reflectivity  of  the  mirrors.  If  the  light  were  to  complete  yet  another 
roundtrip, 


I(aL)  =  I0T2R4e~5aL 


(A.4) 
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The  pattern  now  becomes  distinguishable  enough  to  rewrite  the  expression  as  a  sum¬ 
mation, 


I(aL)  =  I0T2e~a  L  J2(Re~a  L)2n  (A.5) 

n 

If  the  number  of  passes,  n,  is  many  orders  of  magnitude  greater  than  the  length  of 
the  cavity,  L,  then  equation  (A. 5)  can  be  approximated  in  the  infinite  limit  by  an 
integral. 

/»oo 

I(aL)  =  I0T2e~a  L  /  ( Re~aL)2n  dn  (A.6) 

Jo 

One  way  to  reach  the  gain  factor  G  is  to  continue  this  development  by  solving  for  the 
above  integral.  Although  that  would  be  the  simpler  way  to  reach  G,  the  form  would 
be  too  complicated  to  understand  physically  (see  equation  (A. 7)). 


G  = 


2  ln(e 


-aL\ 


(A.7) 


Therefore,  what  needs  to  be  done  is  to  continue  down  a  more  complicated  and  rigorous 
road  in  order  to  find  the  simpler  and  easy-to- visualize  form  of  G  1 .  This  is  done  by 
solving  directly  for  the  summation  in  (A. 5).  Solving  gives, 


-,—cx.L 


I(aL)  =  I0T2 


(1  —  Re 


—2aL\ 


(A.8) 


Which  can  be  re-written  as, 


-,—aL 


I(oiL)  =  I0T2 


(1  -  Re~aL)(  1  +  Re~aL) 


(A.9) 


Now  we  need  to  include  another  intensity  which  is  important  to  multi-pass  absorption 
spectroscopy,  which  is  the  output  intensity  without  the  absorbing  species  present 

1These  series  of  steps  was  carried  out  by  Dr.  Skip  Williams 
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[I(aL  =  0)].  The  next  step  shows  the  beauty  of  obtaining  the  absorption  spectra 
through  direct  attenuation  methods.  If  the  following  relation  is  written, 


I(aL  =  0) 
I(aL) 


(A. 10) 


then  this  gives  the  direct  measurement  of  the  absorbing  signal  is  achieved.  Consider 
for  the  moment  some  fictitious  numbers.  Let  I(aL  =  0)  =  5  (the  intensity  of  the 
output  intensity  without  the  absorbing  species)  and  I{aL )  =  4.9999  (the  output 
intensity  with  the  very  weak  absorbing  species).  If  these  numbers  are  applied  to 
(A.  10),  then  the  result  is  the  intensity  of  the  absorbing  feature,  in  this  case  2  x  10~5. 
Now,  if  equation  (A. 9)  is  applied  to  equation  (A. 10),  then  we  have, 


I(aL  =  0) 
I(aL ) 


Ml 

(l-R)(l+R) 

IoT2e~aL 

( l-Re-aL)(l+Re-aL ) 


(A.ll) 


I(aL  —  0)  _  (1  -Re~aL)(l  +  Re~aL) 

I(aL)  “  (1-  R)(l  +  R)e~aL 


Using  the  following  relations, 


(1  —  J?exp(— aL)) 


(1-R) 


1  + 


R{  1  -  e~aL ) 

(!  -R) 


(A.12) 


(A.13) 


(1  +  Rexp(-aL)) 


(1  +  R) 


R(l-e~aLy 
(1  ~R)  . 


Now  substitute  (A.13)  and  (A.  14)  into  (A.12)  to  obtain, 


I(aL  =  0) 
I(aL) 


1  + 


R{  l-e~aL) 

(1  -R) 


R{  l-e~aL) 

(1  -R) 


o  —  OlL 


(A. 14) 


(A. 15) 


A  term  which  is  known  as  the  single-pass  absorption  can  now  be  defined  as  A  = 
(1  —  e~aL).  This  description  is  more  readily  apparent  when  a  series  expansion  is 
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performed  on  the  exponential  yielding  a  first  approximation,  A  aL.  Therefore, 


I(aL  =  0) 
I (aL) 


1  + 


R) 

(i -R) 


A 


1  - 


R 

(i -R) 


A 


0  —  OlL 


A  -  1 


(A. 16) 


The  next  three  steps  are  simply  algebraic  manipulation. 


I(aL  =  0) 
I(aL) 


1  + 


R 


(1-R) 


A 


R 


(1-R) 


A 


R2 

1-R2 


A2 


D—aL 


(A.17) 


I(aL  =  0) 
I(aL) 


o—a.L 


- 1 


(A.18) 


I(aL  =  0) 
I(aL ) 


1  +  2A  (1^2)  -  (1^2)  A2 
e~aL 


(A.19) 


The  following  step  combines  the  quantity  1  on  the  right-hand  side  by  using  the  com¬ 
mon  denominator.  This  is  done  by  using  the  relation  1  =  e_2L  ■ 


I(aL  —  0)  1  _  (X  _  e  aL)  +  2A  (l5^)  _  (A)  A2 

I(aL)  e~aL 


(A. 20) 


After  combining  the  1,  another  A  =  (1  —  e  aL )  term  appears.  Making  this  substitution 
and  extracting  the  common  eaL  from  the  right-hand  side  gives, 


=  0)  _  t  _  olL 

I{clL ) 


A  +  2A 


R2 


1  -  R 2 


-  A 


R 2 


1-R2 


(A.21) 


The  cavity  enhancement  factor  is  now  defined  as, 


G  = 


2  R2 


1-R2 


(A. 22) 
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Obviously,  equation  A. 22  is  a  lot  more  physically  comprehendible  than  equation  A. 7 
above.  Equation  A.  22  shows  explicitly  how  important  the  reflectivity  of  the  mirrors 
are.  The  closer  R  approaches  100%,  the  larger  G  becomes.  Recall  that  this  experiment 
uses  mirrors  with  reflectivities  of  greater  than  99.997%.  (An  average  value  for  G  in 
this  experiment  was  80,000.  This  yields  reflectivities  of  R  =  99.9988%.)  Substituting 
this  value  of  G  gives, 


I(aL  =  0) 
I(otL) 


A  +  AG  - 


A2G~ 

~Y~ 


Multiplying  through  by  eaL  and  pulling  out  a  factor  of  GA  yields, 


(A. 23) 


I(aL  =  0) 
I(aL) 


GA 


+  eaL 


eaLA 

2 


GA 


aL 


G 


+  e' 


aL 


eaL  (1  -  e~aL) 


(A. 24) 


Now  if  we  expand  out 


eaL(l-e-aL) 

2 


and  combine  terms,  we  end  up  with 


I(aL  =  0) 
I(aL) 


GA 


(A. 25) 


The  next  step  in  this  development  brings  light  as  to  why  this  particular  technique 
(multi-pass  absorption  spectroscopy)  really  only  works  for  very  weak  species.  Since 
the  absorption  cross  section  is  so  small  (  10~24)  then  eaL  — >  1.  Also,  because  we  know 
that  G  will  be  high  (  80,000)  then  ^  — >  0.  Therefore  we  are  left  with  the  following. 


I(aL  =  0) 
I(aL) 


GA 


(A. 26) 


This  is  actually  the  form  used  inside  the  macro  file  of  the  program  used  for  the  data 
analysis  (IGOR  Pro).  This  step  was  already  explained  in  the  Procedure  section  in 
Chapter  III. 

Equation  (A. 26)  is  still  not  the  common  form  used  in  ICOS  theory.  What  is 
now  required  is  a  bit  of  mathematical  messaging.  Multiplying  through  by  I{otL)  in 


equation  (A. 26)  gives, 


I(aL  =  0)  -  I(aL )  =  GA  I(aL) 


(A. 27) 


Next,  we  need  to  add  the  quantity  GA  I(aL  =  0)  to  both  sides. 


I(aL  =  0)  -  I(aL )  +  GA  I(aL  =  0)  =  GA  I(aL)  +  GA  I(aL  =  0)  (A.28) 


Now,  subtract  both  sides  by  GA  I(aL), 


I(aL  =  0)  -  I(aL)  +  GA  I(aL  =  0)  -  GA  I(aL )  =  GA  I(aL  =  0)  (A.29) 


Simplifying  the  left-hand  side, 


(1  +  GA)[I(aL  =  0)  -  I{otL)\  =  GA  I(aL  =  0)  (A.30) 


Next  rearrange  the  terms,  so  that  the  two  GA  terms  are  on  the  right-hand  side. 


I(aL  =  0)  -  I(aL)  GA 

I(aL  =  0)  I  -  GA  {  1 

If  we  define  the  change  in  the  steady-state  cavity  output  as  A /  =  I(aL  =  0)  —  I(aL), 
then  we  have  the  following  usable  form  that  is  the  basis  for  ICOS  theory. 


A  /  _  GA 

I(aL  =  0)  ““  1  +  GA 


(A. 32) 


After  further  analysis  of  this  relation,  it  becomes  clear  that  the  multi-pass 
method  increases  the  intensity  of  the  original  absorption  output  signal  by  a  significant 
amount.  As  mentioned  in  Section  3.1.2,  equation  A. 32  says  that  the  absorption  infor¬ 
mation  (A)  is  contained  within  the  change  in  steady-state  output  of  the  cavity.  This 
allows  for  a  direct  measurement  of  the  absorption  feature  [16].  For  weak  absorption 
(GA  1),  the  cavity  provides  a  linear  absorption  signal  gain.  Physically,  G  refers  to 
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the  number  of  optical  passes  that  occur  within  the  cavity  decay  time  (G  ~  z^)  and 
is  also  simply  related  to  the  cavity  finesse  (G  ~  [12]. 
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Appendix  B.  Term  and  Band  Energies 


Matlab® 


The  following  code  calculates  the  term  energies  for  the  b X3£~  (1,0)  and 
the  a1  Eg  <—  X3£“  (0,0)  transitions.  The  b 1 E+  a1  Ag  (1,0)  band  is  then  found 
by  subtracting  the  appropriate  term  energies.  Lists  for  J  =  1  —  20  are  generated  for 
each  of  these  scenarios.  The  Hamiltonian  and  spectroscopic  constants  used  in  this 
experiment  are  those  found  by  Cheah  and  Lee  [27]. 


B.  1  Matlab  Code 

Listing  B.l:  The  following  is  the  Matlab®  file  used  in  the  thesis. 

( Appendix4  /  cheahter  mband2 .  m) 

1  function  [] =  cheahtermband 

'/.This  program  will  generate  term  energies  for 

y«a(v  =  0)  and  b(v  =  l)  .  The  line  positions  for  the  five  branches 
‘/.are  then  calculated  for  the  b-a  (0,1)  band. 


6 

°/« Spectroscopic  constants  according  to  Cheah  and  Lee  2000 
v0a=7888 . 107842 ; 

B0a  =  1 . 41780381 ; 

D0a=5 . 1 1 144e -6 ; 

11  Del =2 ; 

v0b=13120 .918775 ; 

BOb  =  1 . 39124421 ; 

DOb  =5 . 3641 e -6 ; 

16  Sig=0 ; 

vlb  =  14525 . 65553 ; 

B lb  =  1 . 3729659 ; 

D lb  =5 . 4086  e -6 ; 

21 

f pr intf ( 1 \nTerm  energies  for  02(a),  v=0 ,  in  wavenumber s . \n ’) ; 

‘/.Generates  the  energies  in  wavenumbers  for  02(a),  v=0 
f  or  Ja  =  2 : 20 

26  Fa(Ja)=v0a  +  B0a*(Ja*(Ja  +  1)  -  Del~2)  -  D0a*(Ja*(Ja  +1)  - 

Del ~2) ~ 2 ; 

f  print  f  (’(  J  =°/0 1 . 0  f  )  =  ‘/,9.9e  \n  ’  ,  Ja  ,  Fa  (  Ja)  )  ; 

end 

f pr intf ( 1 \nTerm  energies  for  02(b),  v  =  l  ,  in  wavenumbers ,\n’) ; 

31 

‘/.Generates  the  energies  in  wavenumbers  for  02(b),  v  =  l 
for  Jb  =2 : 2 : 20 
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Fbl(Jb)=vlb  +  B  lb  *  (  Jb  *  (  Jb  +  1)  -  Sig~2)  -  Dlb*(Jb*(Jb  +1)  -  . 

Sig~2)-2; 

f  print  f  (  1  (  J  =°/o  1 . 0  f  )  =  7.9. 9e  \n ’  , Jb , Fb 1 ( Jb  )  )  ; 

36  end 

f pr intf ( 1 \nThe  following  line  positions  are  for  the  b-a  (1,0)  band 
in  wavenumber s \n 1 ) ; 

for  Ja=2 : 20 
41  for  Jb  =  2 : 2 : 20 

if  Jb- Ja==-2 

0  (  Ja) =Fbl ( Jb) -Fa( Ja)  ; 
f  printf  (  ’  0  C/,1  .  Of  )  =  7.9. 9e 

end 

46  if  Jb- Ja==-1 

P( Ja)=Fbl ( Jb) -Fa( Ja)  ; 
fprintf  (  ’P  (7,1  .  Of  )=  %9. 9e 

end 

if  Jb- Ja==0 

51  Q ( Ja) =Fbl ( Jb) -Fa( Ja)  ; 

f  printf  (  ’  Q  (7,1  .  Of  )  =  7, 9 . 9 e 

end 

if  Jb - Ja==l 

R(  Ja)=Fbl ( Jb) -Fa( Ja)  ; 

56  f  printf  (  ’  R  (7.1  ■  Of  )  =  709.9e 

end 

if  Jb- Ja==2 

S( Ja)=Fbl ( Jb) -Fa( Ja) ; 
fprintf  (  ’  S  (  %  1 . 0f)=  7,9. 9e 

61  end 

end 

end 

B.2  The  Theoretical  Term  Energies 

Term  energies  for  O 2(a),  v=0  in  wavenumbers. 

(J=2)=7.890943429e+003 
(J=3)=  7.899449945e+003 
(J=4)=7.910791394e+003 
(J=5)=  7.924967286e+003 
(J=6)=  7.941977006e+003 
(J=7)=  7.961819819e+003 
(J=8)=  7.984494866e+003 
(J=9)=  8.010001165e+003 


\n ’ , Ja , 0 ( Ja) ) ; 


\n ’ , Ja,P(Ja)) ; 


\n ’ , Ja , q ( Ja) ) ; 


\n’ , Ja , R ( Ja) ) ; 


\n ’ , Ja , S ( Ja) ) ; 
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(J=10)=  8.038337614e+003 
(J=H)=  8.069502984e+003 
(J=12)=  8.103495926e+003 
(J=13)=  8.140314969e+003 
(J=14)=  8.179958518e+003 
(J=15)=  8.222424854e+003 
(J=16)=  8.267712139e+003 
(J=17)=  8.315818409e+003 
(J=18)=  8.366741578e+003 
(J=19)=  8.420479440e+003 
(J=20)=  8.477029662e+003 

Term  energies  for  02(b),  v=l  in  wavenumbers. 
(J=2)=1.453389313e+004 
(J=4)=  1.45531 1268e+004 
(J=6)=  1.458331056e+004 
(J=8)=  1.462448104e+004 
(J=10)=  1.467661633e+004 
(J=12)=  1.473970659e+004 
(J=14)=  1.481373985e+004 
(J=16)=  1.489870210e+004 
(J=18)=  1.499457726e+004 
(J=20)=  1.510134713e+004 

B.3  The  Predicted  Band  Energies 

The  following  line  positions  are  for  the  b-a(0,l)  band  in  wavenumbers. 
Q(2)=  6.642949702e+003 
S(2)=  6.662169255e+003 
P(3)=  6.634443185e+003 
R(3)=  6.653662739e+003 
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0(4)=  6.623101736e+003 
Q(4)=  6.642321290e+003 
S(4)=  6.672519163e+003 
P(5)=  6.628145399e+003 
R(5)=  6.658343271e+003 
0(6)=  6.611135679e+003 
Q(6)=  6.641333551e+003 
S(6)=  6.682504031e+003 
P(7)=  6.621490738e+003 
R(7)=  6.662661218e+003 
0(8)=  6.598815691e+003 
Q(8)=  6.639986171e+003 
S(8)=  6.692121469e+003 
P(9)=  6.614479871e+003 
R(9)=  6.666615169e+003 
0(10)=  6.586143423e+003 
Q(10)=  6.638278721e+003 
S(10)=  6.701368973e+003 
P(ll)=  6.607113351e+003 
R(ll)=  6.670203603e+003 
0(12)=  6.573120409e+003 
Q(12)=  6.636210660e+003 
S(12)=  6.710243923e+003 
P(13)=  6.599391617e+003 
R(13)=  6.673424880e+003 
0(14)=  6.559748069e+003 
Q(14)=  6.633781332e+003 
S(14)=  6.718743587e+003 
P(15)=  6.591314995e+003 
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R(15)=  6.676277251e+003 
0(16)=  6.546027711e+003 
Q(16)=  6.630989966e+003 
S(16)=  6.726865117e+003 
P(17)=  6.582883696e+003 
R(17)=  6.678758847e+003 
0(18)=  6.531960527e+003 
Q(18)=  6.627835678e+003 
S(18)=  6.734605553e+003 
P(19)=  6.574097817e+003 
R(19)=  6.680867691e+003 
0(20)=  6.517547595e+003 
Q(20)=  6.624317469e+003 
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Appendix  C.  Analyzed  Data 

This  appendix  contains  all  of  the  21  lines  observed  in  this  research  effort.  Please 
reference  the  left-axis  (which  describes  the  signal  intensity)  for  each  line,  because 
these  lines  are  not  scaled  relative  to  each  other.  Instead,  each  line  has  been  zoomed 
in  to  fill  the  graph.  These  particular  graphs  were  the  one’s  used  to  measure  the 
absolute  line  positions  to  within  0.004cm~1.  The  helium  and  oxygen  flow  rates  were 
4000  seem  and  400  seem  respectively,  and  the  average  pressure  within  the  optical 
cavity  was  23  torr. 
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Figure  C.l:  The  0(6)  line  position  at  6611.1418cm  1. 
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Figure  C.2:  The  0(8)  line  position  at  6598.8226cm  1. 


Figure  C.3:  The  P(9)  line  position  at  6614.4826cm  1. 
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Figure  C.4:  The  P(13)  line  position  at  6599.3977cm  T 
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Figure  C.5:  The  P(15)  line  position  at  6591.3228cm  1. 
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Figure  C.6:  The  Q(2)  line  position  at  6642.9513cm  1. 
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Figure  C.7:  The  Q(4)  line  position  at  6642.3227cm  1. 
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Figure  C.8:  The  Q(6)  line  position  at  6641.3324cm  1. 
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Figure  C.10:  The  Q(10)  line  position  at  6638.2786cm-1. 


Figure  C.ll:  The  Q(12)  line  position  at  6636.2118cm  T 
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Figure  C.12:  The  Q(14)  line  position  at  6633.7802cm  1. 


Figure  C.13:  The  Q(16)  line  position  at  6630.9987cm  1. 
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Figure  C.14:  The  R(3)  line  position  at  6653.6667cm  F 
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Figure  C.15:  The  R(7)  line  position  at  6662.6576cm  F 
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Figure  C.16:  The  R(9)  line  position  at  6666.6174cm  F 


Figure  C.17:  The  R(ll)  line  position  at  6670.2048cm 
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Figure  C.18:  The  R(13)  line  position  at  6673.4308cm  1. 


Figure  C.19:  The  R(17)  line  position  at  6678.7650cm  F 
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Figure  C.20:  The  R,(19)  line  position  at  6680.8693cm  1. 
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Figure  C.21:  The  S(8)  line  position  at  6692.1316cm  1 . 
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Cross  Section  (cm ) 


Appendix  D.  Developing  the  Hamiltonian 


This  section  describes  the  various  approximations  used  in  quantum  mechanics 
which  lead  up  to  the  rotation-electronic  interaction  approximation  listed  in 
Hund’s  cases  (a-e).  The  Hamiltonian  for  Hunds  case  (a)  will  then  be  briefly  mentioned, 
since  this  was  the  Hamiltonian  used  by  Cheah  and  Lee  [27]  to  come  up  with  their 
term  energies  b 1S+  <—  and  a1Ag  <—  A"3E“,  and  ultimately  was  used  in  this 

experiment  to  derive  the  predicted  line  positions  for  b <—  a1Ag.  Recall  that 
the  case  (b)  Hamiltonian  describes  the  blY,^  <—  a1Ag  transition  directly,  however 
spectroscopic  constants  do  not  yet  exist,  so  the  case  (b)  could  not  be  used  in  this 
experiment  to  determine  the  line  positions. 


D.l  The  Rigid  Rotor  to  the  Case  (a)  Hamiltonian 

Considering  rotation  only,  a  good  first  order  approximation  of  the  eigenenergies 
produced  by  the  homonuclear  diatomic  is  to  treat  it  as  a  rigid  rotor.  The  molecule 
can  then  be  seen  as  two  like  masses  connected  by  a  rigid  rod.  When  rotation  occurs 
the  eigenenergies  are  simply  [8], 


h2J(J  +  1) 
87 Agrl 


(D.l) 


Where  g  is  the  reduced  mass  of  the  molecule  and  re  is  the  distance  between  the  two 
atoms.  If  the  term  energies  are  defined  as  F(J)  =  ^  and  the  moment  of  inertia 
defined  as  I  =  gr2,  then  the  rotational  constant,  B,  can  be  written  as 


e  87 r2c/ 


(D.2) 


However,  molecules  are  not  rigid,  so  a  better  model  must  now  be  considered. 
If  instead  the  atoms  of  the  diatomic  are  connected  by  a  spring,  then  this  describes 
the  non-rigid  rotor  approximation.  In  this  configuration,  any  increase  in  rotation  will 
achieve  an  increase  in  it’s  moment  of  inertia.  As  a  result,  centrifugal  distortion  occurs. 
The  centrifugal  distortion  constant  De  is  dependent  on  the  equilibrium  vibration 
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frequency,  uje  [8]. 


4  Bl 


(D.3) 


Equation  D.l  then  becomes, 


F{J)  —  BeJ(.J  +  1)  —  De[J(J  +  l)]2 


(D.4) 


Equation  D.4  adequately  describes  the  energy  levels  of  a  rotating  diatomic  whose 
potential  energy  curve  is  modeled  by  the  harmonic  oscillator.  However,  when  vibration 
is  present,  the  distance  between  the  two  atoms  in  the  diatomic  increase  and  decrease 
from  their  equilibrium  position.  The  harmonic  oscillator  is  only  a  good  approximation 
for  energy  levels  at  or  close  to  the  equilibrium  position.  Real  molecules  must  therefore 
be  modeled  after  the  anharmonic  oscillator. 

Considering  vibration  only,  the  eigenvalues  for  the  anharmonic  oscillator  are  [8], 


hcOJe (V  +  -)  -  hcQJeXe(v  +  W  +  hajeye(v  +  -)3 


(D.5) 


Where  uJeXe  and  a )eye  are  simply  coefficients  to  the  quadratic  and  cubic  terms,  and  v 
is  the  vibrational  state.  Writing  equation  D.5  in  terms  of  term  energies  (G(v)), 


G(v)  =  u>e  (  v  + 


1 

2 


a^eXe  T  ^ eVe 


(D.6) 


Another  term  which  is  used  often  in  spectroscopy  is  vq. 


Vo  =  G(y')  -  G(v") 


(D.7) 


This  is  the  vibrational  term  that  describes  the  energy  difference  between  an  upper 
and  lower  vibrational  state. 


An  even  more  realistic  model  of  the  diatomic  is  the  vibrating  rotator,  which 
involves  vibration,  rotation,  and  the  interaction  between  them.  If  there  didn’t  exist 
an  interaction,  the  corresponding  energy  of  the  vibrating  rotator  would  simply  be  the 
sum  of  equations  D.4  and  D.6.  However,  during  vibration  the  internuclear  distance 


will  change,  which  then  changes  the  molecule’s  moment  of  inertia,  and  in  turn  changes 


the  molecule’s  rotational  constant  Be  of  equation  D.2.  Since  the  period  of  vibration 
is  very  small  compared  to  that  of  the  period  of  rotation,  a  mean  Be  value  for  the 
rotational  constant  in  the  vibrational  state  can  be  used.  In  other  words,  as  v  increases 
the  molecule  spends  more  of  it’s  time  at  large  r  where  the  potential  energy  curve 
is  flatter.  Thus  the  average  internuclear  separation  (r)  increases  with  v,  while  Bv 
decreases  [8]. 


(D.8) 


The  above  equation  is  parameterized  using  wave-mechanical  calculation  to  obtain  [8], 


(D.9) 


Where  ae  and  ye  are  just  coefficients.  The  same  process  is  performed  for  the  centrifu¬ 
gal  distortion  constant, 


(D.10) 


Where  /3e  is  also  a  coefficient.  Substituting  the  constants  at  equilibrium  for  their 
vibration- rotation  interaction  counterparts,  equation  D.4  can  then  be  expressed  as, 


F(J)  =  BVJ ( J  +  1)  -  DV[J(J  +  l)]2 


(D.ll) 


There  is  still  another  model  which  more  accurately  describes  the  energy  levels 
for  a  homonuclear  diatomic,  because  the  interaction  between  the  rotational  and  elec- 


tronic  motions  must  also  be  conidered.  Hund  was  the  first  to  categorized  the  different 
coupling  schemes  between  the  rotational  and  electronic  motions.  In  fact,  Hund  dis¬ 
covered  that  there  are  five  ways  for  rotational  and  electronic  motions  to  couple  to 
each  other.  These  are  the  five  well  known  cases  (a-e). 

Hund’s  case  (a)  states  that  the  electronic  motion  of  the  molecule,  which  includes 
both  the  orbital  and  spin  angular  momentum,  is  strongly  coupled  to  the  internuclear 
axis  (the  axis  of  rotational  symmetry).  This  coupling  scheme  is  very  similar  to  the 
symmetric  top.  Therefore  the  ro-vibronic  energy  of  a  molecule  in  Hund’s  case  (a) 
has  the  form  of  the  symmetric  top.  Returning  to  equation  D.l  and  including  the 
symmetric  top  energy, 


F(J)  =  BeJ(J +1)  +  (A- Be)A2  (D.12) 

Where  HA2  can  be  understood  as  the  spin-orbit  coupling  term  [25],  but  since  S=0  for 
singlet  states  (i.e.  and  XA),  then  A=0.  Therefore,  equation  D.12  reduces  to, 

F(J)  =  BV[J(J  +  1)  —  A2]  (D.13) 

The  final  Hamiltonian  used  by  Cheah  and  Lee,  includes  the  symmetric  top  (equa¬ 
tion  D.13,  the  vibrating  rotator  D.ll  and  the  vibrational  term  D.7  [27], 

F(J)  =  u0  +  B0[J(J  +  1)  -  A2]  -  D0[J{ J  +  1)  -  A2]2  (D.14) 
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14.  ABSTRACT 

Off-axis  integrated-cavity-output  spectroscopy  (ICOS)  has  been  applied  to  detect  singlet  oxygen  in  the  (1,0)  band  of  the 
electric-quadrupole  allowed  Noxon  system  < —  a1  A.g)  of  oxygen.  Absolute  line  positions  were  determined  for  21 

lines  to  within  0.004C77W1 .  Ten  of  these  lines  were  pressure  broadened  from  25  to  100  torr,  and  pressure  broadening 
coefficients  were  determined  from  =  3.61  X  10_6C7Tt_1  ftorv  to  5.81  X  10 cmT1  / torr .  The  singlet  oxygen  was 
generated  in  a  microwave  plasma  and  the  off-axis  ICOS  measurement  system  consisted  of  an  82-cm  long  high-finesse 
(~  248000)  optical  cavity.  The  mirror  reflectivity,  R,  was  determined  using  a  cavity  ringdown  measurement  and 
observing  ringdown  times  of  220  /is  in  a  range  from  1494  nm  to  1512  nm.  The  absorption  within  the  cavity  leads  to  a 
decrease  in  transmitted  intensity,  A/,  given  by  A I / IQ  =  GAj  (1  +  GA ),  where  IQ  is  incident  laser  intensity,  A  is  the 
single  pass  absorption  (~  5  X  1CD8),  and  G  =  2R2/(1  —  R2)  ~  80,  000  is  the  cavity  enhancement  factor. 

15.  SUBJECT  TERMS 

cavity  enhancement  factor,  cavity  ringdown,  electric-quadrupole  transition,  off-axis  ICOS,  singlet  delta  oxygen 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF 

PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

Glen  P  Perram 

u 

u 

u 

uu 

111 

19b.  TELEPHONE  NUMBER  (include  area  code) 

(937)  255-3636,  ext  4504 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


